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(Received April 29, 1946) 


| This paper presents data on the internal friction of four 
single crystals of zinc while oscillating longitudinally, and 
a description of various slow speed tension tests on a fifth 
crystal within and beyond the elastic limit. In the first 
named measurements the behavior of the crystals bears 
little resemblance to that of crystals of zinc of greater 
purity prepared by another method, which has previously 
been reported (see reference 1). The most outstanding 
feature is that the decrement, although higher at the lowest 
stress amplitude than for the previous ones, shows very 
little rise with increasing stress amplitude, even up to 
stresses far beyond the statically-determined elastic limit. 
The difference in behavior seems to be caused by the 
difference in purity of the zincs. An optically mosaic 


structure, such as that described by Schilling (see refer- 
ence 4), does not appear to be responsible. In the slow 
speed tension tests the zinc crystal, after a period of self- 
annealing at room temperature, has a Hooke’s law region 
up to about 70 g/mm? (R.S.S.=28 g/mm?) and at greater 
stress shows transient and steady creep. Under some 
circumstances the transient creep is in two parts, a slow 
starting creep followed by much more rapid creep. The 
crystal, when measured following a permanent strain, 
shows creep at loads well below the previous elastic limit 
and has a large amount of hysteresis. Self-annealing occurs 
during rest periods and the initial elastic limit is gradually 
regained while the hysteresis disappears. 





_ INTRODUCTION 


HE only previous work on the internal fric- 

tion of zinc single crystals is that described 

by one of the authors in a previous paper.! The 
specimens used in this earlier investigation were 
grown from spectrographically pure zinc (99.999 
percent Zn) by the Bridgman method in Pyrex 
glass tubes coated with graphite. The damping 
behavior of these crystals was very complicated, 
but may be briefly summarized by saying that 


* Westinghouse Research Fellow. Now at the Ord- 
nance Laboratory, Frankford Arsenal, Philadelphia, Penn- 
sylvania. 

1T, A. Read, Phys. Rev. 58, 371 (1940). 


after a one- or two-day anneal at room tempera- 
ture their internal friction was very low (decre- 
ment about 1X10-°) at small amplitudes of 
vibration, but increased with increasing ampli- 
tude. This increase was very rapid for crystals 
whose hexagonal axis made an angle of from 15 
to 75 degrees with the axis of the rod. For ex- 
ample, the decrement of a crystal for which this 
angle was 74.5° was 1.7X10~ at a stress ampli- 
tude of .15 g/mm?, and 144X10- at a stress 
amplitude of 22.2 g/mm?. Both these figures for 
the stress refer to the values at displacement 
nodes. The corresponding values of the strain 
amplitude are 1.210~-* and 176X10-5. 
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The objects of the investigation described in 
the present paper were to compare the internal 
friction of zinc single crystals grown by another 
method with that of the Bridgman crystals, and 
to study the relation between internal friction 
and behavior in slow tension tests. 


SPECIMEN MATERIAL 
The zine crystals used in the present investi- 
gation were grown from ‘‘Bunker Hill’ zinc, 
a commercial zinc whose total impurity content 


TABLE I. Description of crystals. 





Designation Orientation R.S.S. factor Remarks 
2 45° 0.25 O.M. 
3 a 0.29 G.C. Pb, Sn 
4 70° 0.32 G.C. 
5 81° 0.17 O.M. Pb, Sn 
6 63° 0.41 G.C. 


O.M.: Optically mosaic. 
G.C.: Good crystal. 
Pb, Sn: Lead and tin impurities. 


is about 0.01 percent. It is to be noted that this 
purity is somewhat less than that of the spectro- 
graphically pure zinc used in the earlier investiga- 
tion referred to above. At the conclusion of the 
work a spectrographic analysis was made on four 
of the crystals and on a small ingot cast from 
the original material. It revealed that crystals 
Nos. 2 and 4 were identical with the original 
material but that Nos. 3 and 5 contained excess 
amounts of lead and tin,? which had apparently 
been introduced accidentally by the use of a 
‘“‘seed”’ crystal contaminated with these metals. 
The effect of the difference in purity between 
the two pairs of crystals was very marked as 
may be seen below. 

These crystals were grown by the modified 
Kapitza method described by Cinnamon.’ The 
molds were made of well-baked transite with 
trapezoidal grooves, 0.25 cm? in cross section 
for crystals Nos. 2-5 and 0.5 cm? for crystal 
No. 6. The surface of specimens grown this way 
is somewhat rougher than that of crystals grown 
in glass molds. This is caused by the wrinkling 
of a very thin surface film of oxide when the 


* In crystals Nos. 2 and 4 and the ingot Sn was either 
not detected or less than 0.001 percent, while Pb was 
either not detected or less than 0.002 percent. In Nos. 3 
and 5 Pb was about 0.01 percent, and Sn 0.1 percent 
(No. 3) and 0.03 percent (No. 5). 

?C. A. Cinnamon, Rev. Sci. Inst. 5, 187 (1934). 
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metal contracts in freezing. The method is 
believed, however, to be more likely to produce 
more nearly perfect crystals because of lack of 
constraint on the growing crystal. Two of the 
crystals were, however, purposely grown of the 
“optical mosaic’ type.*® 

The specimens are described in Table I. The 
orientation is given as the angle between the 
hexagonal axis and the specimen length. The 
third column gives the factor by which a longi- 
tudinal stress must be multiplied to obtain the 
resolved shear stress in the basal (glide) plane. 
Further resolution into the glide direction would 
require an additional factor which must, how- 
ever, always lie between 1.0 and 0.866. Any 
deviation of this from 1.0 has been neglected in 
the following. 

The decision as to the presence of an optically 
mosaic structure was made on the basjs of 
examination of the portions of the crystal ad- 
jacent to the specimen, and of the specimen 
itself after the completion of the measurements. 

After growth the first four crystals were cut 
to approximately the proper length for one half- 
wave of longitudinal vibration with a jeweler’s 
hack saw run automatically and so arranged 
that the blade pressed very lightly against the 
crystal. They were then packed in Vaseline in 
brass tubes which were in turn packed in cotton 
in a sturdy box. This box was mailed from Iowa 
City to East Pittsburgh. The authors feel that 
the care taken in this handling of the specimens 
sufficed to insure that their internal friction was 
not affected thereby. The results of the measure- 
ment are consistent with this view with one 
exception, as noted below. In East Pittsburgh 
the crystals were cut to exact length with a 
jeweler’s saw after being waxed in two V-blocks 


‘ with beeswax and rosin. One end of each crystal 


was ground flat by holding the crystal lightly 
in a V-block and grinding it on a cast iron surface 
plate. This same method of preparation was em- 
ployed in the previous work.' Crystal No. 6 was 
mounted in the measuring apparatus as grown. 
The crystal was 15 cm long and measurements 
were made on a central 9 cm length. 


4H. K. Schilling, Physics 6, 111 (1935). 

5 Electrical and thermal resistivities of optically mosaic 
crystals differ in an erratic way from those lacking this 
structure. W. J. Poppy, Phys. Rev. 46; 815 (1934) and 
C. A. Cinnamon, Phys. Rev. 46, 215 (1934). 
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INTERNAL FRICTION 


The internal friction of these crystals was 
measured for longitudinal vibration at a fre- 
quency of 33.5 kilocycles per second by the 
composite piezoelectric oscillator method. The 
apparatus used was similar to that previously 
described.! The dependence of the internal fric- 
tion of the specimen crystals on the amplitude 
of vibration was obtained by measuring the 
impedance of the composite piezoelectric oscil- 
lator as a function of the voltage impressed on 
electrodes attached to the quartz and of fre- 
quency in the neighborhood of the resonant 
frequency. It can be readily shown that the stress 
amplitude o at a displacement node in the 
specimen is given by the formula 

2E&(10)'K 
o =———— X 10? dynes/cm?’, (1) 

RoV 





where E= Young’s modulus in dynes/cm? of the 
specimen in the direction of its length. 


&=r.m.s. voltage applied to electrodes. 

K =proportionality constant between equivalent induct- 
ance of composite piezoelectric oscillator and its 
mass in henrys per gram. 

Ro=resistance at resonance in ohms for given value of &. 

\’=velocity, in cm/sec., of sound in the specimen in the 
direction of its length. 


The formula (1) was given in reference 1 in a 
somewhat less convenient form. 
For all of the measurements the specimens 
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Fic. 1. The internal friction of the four crystals investi- 
gated as a function of the stress amplitude at a displace- 
ment node. The arrows indicate a change with time at 
constant driving voltage. The crystal number is indicated 
on each curve. 
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were suspended in a vacuum and maintained 
at a temperature of 23°C. 

A convenient measure of internal friction is 
the decrement, defined as the ratio of the energy 
dissipated in the specimen in a cycle of vibration 
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Fic. 2. The lower curve gives the internal friction of 
crystal number 4 after a one hour anneal at 330°C. The 
upper curve, repeated from Fig. 1, shows the original 
internal friction. 


to twice the total energy stored in the vibration 
of the specimen. 

The dependence of the internal friction of 
each of the four specimens on the amplitude of 
vibration is given in Fig. 1. These data, obtained 
without any annealing or other treatment except 
as described above, are remarkable for five 
reasons. First, the recovery period of a day or 
two after mounting, during which the previously 
measured crystals! exhibited a decrease in in- 
ternal friction by a factor of from 10 to 100, was 
completely absent for these crystals, although 
exactly the same technique of affixing the speci- 
mens to the quartz rod was used. Second, the 
internal friction of the two crystals of lesser 
purity is only slightly greater at stress amplitudes 
of several hundred g/mm?* than at very small 
amplitudes of vibration. Third, there is no 
marked difference between the behavior of the 
crystals of the same purity which have an opti- 
cally mosaic structure and those which do not 
(compare Nos. 2 and 4 and Nos. 3 and 5). 
Fourth, these crystals have considerably higher 
internal friction at small amplitudes of vibration 
than did the zinc crystals described in reference 1. 
Fifth, the change of the internal friction with 
time, at constant driving force, which occurs at 
the largest amplitudes is a decrease. This change 
is indicated by the arrows. 

Data were obtained on three methods of 
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Fic. 3. (a) The internal friction of crystal number 3 as 
a function both of increasing and decreasing amplitude of 
vibration. The sequence of the measurements is indicated 
by the arrows. (b) The changes in the natural frequency 
of longitudinal vibration of crystal specimen number 3 
which are observed simultaneously with the decrement 
values plotted in Fig. 3(a). 


changing the internal friction of these crystals, 
viz., etching, annealing, and cold-working. 
Etching one of these crystals does not appreci- 
ably affect the dependence of internal friction on 
amplitude, but does produce a lowering of the 
curve as a whole. For example, a three-minute 
immersion of crystal No. 3 in 30 percent HCl 
decreased the small-amplitude decrement from 
7.6X10-° to 4.0 10-5. This treatment decreased 
the weight of this crystal by 0.05 percent. 
Subsequent electrolytic etching of this crystal 
lowered the small-amplitude decrement further 
to 2.3X10~-°. The corresponding loss in weight 
was 3.8 percent. It may be concluded that the 
greater part of the original internal friction of 
this crystal at small amplitudes was a surface 
effect rather than a property of the material. 
As shown in Fig. 1 (and in Fig. 2, where a 
more extensive curve is plotted), the internal 
friction of crystal number 4 was considerably 
higher and rose more rapidly with increasing 
amplitude than for the other crystals. For this 
reason this crystal was selected for an investiga- 
tion of the effect of a heat treatment; it was 
suspected that this crystal may have been 
strained by the handling described above. The 
change in the internal friction which was pro- 
duced by an anneal of one hour at 330°C is shown 
in Fig. 2, with the lower curve depicting the 
results obtained after the anneal. Of the four 
crystals described in this paper, this particular 
crystal exhibited the greatest similarity in proper- 
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ties to the crystals reported on in reference 1. 
For example, the internal friction of this crystal 
always increased with time when any change 
was observed at constant amplitude of vibration. 

Cold-working raises the internal friction of 
these zinc crystals, except in certain cases in 
which the cold-work is introduced by the internal 
friction measurement itself. This exception has 
been mentioned in connection with Fig. 1. The 
more typical case, in which the cold-working 
associated with the internal friction measure- 
ment raises the internal friction is exemplified 
by Fig. 3 (a). In this figure is plotted a sequence 
of measurements on crystal No. 3, beginning at 
a very small amplitude of vibration, proceeding 
in steps to a stress amplitude of 420 g/mm?, and 
then decreasing in steps to the initial amplitude 
of vibration. This curve is similar in shape to 
one given in reference 1, except that here the 
internal friction actually rises to a maximum as 
the amplitude of vibration is decreased. A period 
of only twenty minutes elapsed during the 
measurements. After approximately 24 hours the 
lower branch of the curve may be reproduced. 
The internal friction of the crystal shows a 
more rapid rise than that shown in Fig. 1 
because shortly before making the measurements 
for Fig. 3 this crystal was cold-worked statically 
in the manner described below. 

As is described at greater length in reference 1, 
there is in general associated with the change in 
the decrement produced by a change of the 
amplitude of vibration a change in the natural 
frequency of longitudinal vibration of the speci- 
men rod. This change in the resonant frequency 
is, in many cases, related in a simple way to the 
change in internal friction. In the instances cited 
there a plot of successive values of the change in 
natural frequency against the corresponding 
values of the decrement for a series of amplitudes 
of vibration produced a straight line. That this 
simple relation does not exist for crystal number 
3 under the present conditions of measurement 
is demonstrated by a comparison of Fig. 3 (a) 
with Fig. 3 (b) which gives the changes (de- 
creases) in the natural frequency of this crystal 
which were observed simultaneously with the 
decrement changes plotted in Fig. 3 (a). It is, 
however, important to note that a maximum 
slope of the frequency change curve given in 
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Fig. 3 (b) occurs at the same amplitude of 
vibration as the maximum of the internal friction 
curve of Fig. 3 (a). 

The effect on the internal friction of the cold- 
work introduced by a static load is shown by 
the difference between the lower and upper 
curves of Fig. 4, which were taken, respectively, 
immediately before and after subjecting crystal 
No. 3 to a compressive stress of t.8 g/mm? for 
one minute. The greater part of this increase in 
internal friction disappeared in the course of 
two or three days. These data suggest that a 
previously applied static load which acts for an 
appreciable time is more effective in causing the 
decrement to show the rapid rise with increasing 
stress amplitude (i.e., the specimen acts less 
elastically) than is a much greater stress which 
is rapidly alternated. The slow speed tension 
tests described below seem to afford some clue 
to this behavior. 

SLOW SPEED TENSION TESTS 

Various slow speed tension tests were made on 
crystal No. 6 with the apparatus which has 
been described briefly elsewhere.* The tests con- 
sisted in observing the elongation under load for 
various conditions of the crystal, which were 
determined by the nature of the preceding runs, 
periods of rest between runs, etc. Creep curves 
under steady loads of various values were also 
obtained. Representative experiments are de- 
scribed in detail below. 

The first run taken on the crystal immediately 
after mounting it in the apparatus showed slight 
hysteresis and a lower elastic limit’ than was 
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Fic. 4. The increase in the internal friction of specimen 
number 3 produced by a static compressive stress of 
1.8 g/mm* applied for one minute. Note that this static 
stress is much smaller than the stress amplitude in the 
internal friction measurements. 


6]. H. Swift and E. P. T. Tyndall, Phys. Rev. 61, 359 
(1942). 
7 Ata stress of 43 g/mm’. 
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STRAIN x 10° 


Fic. 5. Crystal No. 6: curve I, in the elastic state. The 
last load produced a permanent set of Al/l=2X10~ 
(not shown on curve). Curves II, III, crystal in inelastic 
state. Curve IV, after two days rest. Curve V, hysteresis 
loop produced by small plastic strain at end of curve IV. 


expected. A slight permanent set (Al//=1.5 
<10-*) was given the crystal at this time. It 
was therefore allowed to rest under the minimum 
load (equivalent to a stress of about 7 g/mm?*) 
for two weeks. It was then found to be in an 
“elastic state’ (see later) and the Young’s 
modulus was carefully measured, using the 
photoelectric multiplying device® and, for safety, 
with the maximum load kept below the pre- 
viously found elastic limit. After a total period 
of 42 days in which the crystal had either rested 
or been subjected to only small loads, the load 
was gradually increased. The data obtained are 
shown in curve I, Fig. 5. It will be seen that 
the proportional region goes to a stress of 67.2 
g/mm?*. The next load (stress=69.9 g/mm?*) 
caused creep to an “‘off-scale”’ deflection and the 
specimen was hurriedly unloaded. The creep, 
however, was observed to occur rather slowly 
at first and then with increasing rapidity. In 
fact the point plotted at 69.9 g/mm* corresponds 
to the scale reading a second or two after com- 
pleting the loading. This phenomenon of slow 
creep followed by much more rapid creep was 
observed many times in subsequent experiments. 
The duration of the slow creep is far too long to 
be accounted for by the inertia of the load or of 
the moving parts of the observing apparatus. 
The permanent set left after unloading was 


® The loading is done by running mercury into a con- 
tainer hung from the lower end of the crystal. About 20 
seconds elapsed between the 67.7 and 69.9 g/mm? readings. 
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Fic. 6. Crystal No. 6: curve VI, crystal in elastic state 
after ten days rest. Permanent set at end of curve VI, 
Al/l=1.7X107* (not shown). Curve VII hysteresis loop 
and creep, immediately following curve VI. Curves VIII 
and IX, one and three days after curve VII. Curve X, 
typical hysteresis loops of crystal in inelastic state. 


Al/l=2X10-. The data shown in curve II 
were taken next.’ This curve shows creep at 
two values of stress (41.4 and 46.9 g/mm*) 
which were previously within the elastic or 
proportional region. The time interval between 
the two points plotted at 41.4 g/mm®? is 3 
minutes; at 46.9, it is 2 minutes. On unloading 
there is a slight recovery in 1 minute as shown 
by the two points at the minimum stress. 
Curve III followed immediately, a well-defined 
and definite hysteresis loop. While at each load 
one point only is plotted, it should be mentioned 
that, if given time, noticeable creep will occur at 
the higher stress points on the increasing stress 
branch. On the descending branch, however, 
no observable creep occurs. After curve II] was 
completed, the crystal was loaded to 50 g/mm* 
and various observations of creep made. It was 
then left under the minimum load for 22 hours, 
following which curve IV resulted. The first, 
straight part of the curve is, within the accuracy 
of observation, identical in slope to the “‘elastic’”’ 
curve (1). Definite creep occurs, however, at the 
highest stress shown (2 minutes interval between 
the two readings), and the following run, curve V, 
again shows hysteresis. This hysteresis is not as 

® The first point on curve II and on the other curves is 
plotted at the correct stress value but the strain is set 
arbitrarily in order to separate or compare certain curves. 
Thus II and III are set together with the same initial 


strain whereas in fact III started from the last point ob- 
served in II. Curves IV and V are similarly arranged. 
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much as it is in curve III and the average slope 
is nearer that of the elastic curve. 

The results depicted in Fig. 6 are very similar 
to those just described. Curve VI is an elastic 
curve taken after the crystal had had quite a 
few runs taken on it but had then rested for 10 
days. Points were taken on this curve for both 
increasing and decreasing loads below 60 g/mm*, 
and no observable hysteresis was found. At the 
highest stress shown the results described in 
connection with curve | were duplicated. The 
permanent set after urloading was Al/]=1.7 
x 10-4. The lower part (hysteresis loop) of curve 
Vil was taken next, and then the complete 
curve with repetition of the upper branch. In 
this curve the time intervals between successive 
points at the same stress are 1 minute, except 
for 2 minutes at 40.2 g/mm/?. After taking the 
upper part of curve VII and returning to the 
minimum load, the lower part (hysteresis loop) 
was duplicated. Curves VIII and IX were taken 
one and three days, respectively, after the data 
for curve VII, with no other runs intervening. 
The gradual recovery of the crystal is obvious. 
In both these curves the ascending branch is 
almost straight and parallel to the ‘‘elastic”’ 
curve (VI). Curve X shows three successive 
runs on a hysteresis loop like that of curve VII 
with points indicated by circles (O), crosses (+), 
and triangles (A), and taken in that order. 
These observations within the 
experimental accuracy and a definite cyclic 
state seems to be indicated. 

Creep data of the same type as that shown by 
the succession of points at each stress in curve 
VII, are represented in a different way in Fig. 7. 
These data -were taken when the crystal was in 
a strained condition which was produced by 
starting with the crystal in an elastic state and 
then going beyond the elastic limit to a perma- 
nent set of Al/l=0.8X10-*. For each curve, 
except XI, the strain at the point marked A is 
arbitrarily plotted as zero although it is really 
the strain reached at the end of the preceding 
run (curve of next lower number). At A the 
increment to the load commenced and was com- 
pleted several seconds before taking the reading 
at B, which is 20 seconds after A. The pronounced 
“‘toe’”’ on curve XV shows the initial slowness 


sets of agree 
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of the transient creep, to which attention has 
been directed previously. 

From the foregoing one may conclude that a 
zinc crystal of this orientation,'® when sufficiently 
rested (self-annealed), has an elastic limit which 
is very well defined by the simultaneous occur- 
rence of the end of the Hooke’s law region and 
the point at which marked and rapid creep 
occurs. This limit is of the type previously re- 
ported for lead crystals* but the subsequent 
behavior of the zinc is very different from that 
of lead when the limit is passed and creep occurs. 
The zinc is now in an inelastic state in which 
Hooke’s law fails both for increasing and de- 
creasing stress and in which marked hysteresis 
occurs. It must be emphasized that the hysteresis 
loop cannot be wholly accounted for by assuming 
creep while the load on the specimen is being 
increased since one should then find, at the end 
of the cycle, a permanent set equal to the 
accumulated deviation from the elastic curve. 
This inelastic state gradually disappears with the 
self-annealing of the zinc. Before complete re- 
covery the material may return to almost the 
original elastic state, with no or very little 
hysteresis, but with an elastic limit markedly 
lower than in the completely recovered state. 
Something like seven or eight days rest at room 
temperature seem to be needed for the limit to 
rise to the value of curves | and VI. 

During the whole set of experiments on crystal 
No. 6 the total accumulated permanent set was 
only Al//=11X10~', or about 0.1 percent. The 
specimen showed no strain nor slip bands on 
its surface. A careful measurement of Young’s 
modulus after an appropriate rest interval yielded 
fk =11.6X10" dynes per cm*®. This is about 3 
percent lower than the initial ‘value of 12.0 
X10" dynes/cm?, corresponding to the crystal 
in the state of curve I, Fig. 5." 

More recent work at the University of Iowa has con- 
firmed the behavior described above for crystals of orienta- 
tions above 45° but has shown that crystals of orientation 
lower than 45° act somewhat differently. 

The critical reader will note that the slope of this 
curve does not correspond to the stated value of Young’s 
modulus. This is because the curves were all plotted with 
the stress computed from a provisional value of cross 
section. The cross section was only obtained accurately at 
the conclusion of the work by cutting out and weighing 
the actual 9-cm length on which observations were made. 
It did not seem worth while to correct all the curves since 


the error in stress is small and is the same for all the 
curves. 
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Fic. 7. Crystal No. 6: creep under successively greater 
loads. At A the load is that shown on the curve of next 
lower number. Between A and B the load is raised to the 
value shown on the curve and thereafter kept constant. 
The loading is complete several seconds before B. The 
interval AB is twenty seconds. Time increases from left 
to right with the lowest point on each curve at t=0. The 
initial strain for each curve is arbitrarily set equal to zero. 


DISCUSSION 


Some of the hitherto puzzling phenomena ob- 
served in connection with internal friction of 
zinc crystals may receive provisional and qualita- 
tive explanations in terms of the slow speed 
observations. If a zinc crystal may be put into 
the inelastic state by accidental strains incurred 
during handling, mounting for measurement, 
etc., then such a crystal would be expected to 
have a high decrement (as measured at the 
minimum stress amplitude). With some rest this 
decrement should decrease, but it might still 
rise with increasing stress amplitude until rested 
sufficiently for a complete return to the elastic 
state. When this state is reached the decrement 
at low stress amplitudes should be a minimum 
and the decrement curve should stay substan- 
tially level until at least the static elastic limit 
is reached. Indeed this point might be exceeded 
somewhat since it is the maximum stress ampli- 
tude which is plotted and this exists in only a 
very small portion of the crystal with less stress 
elsewhere: It is, however, rather obvious that in 
these experiments the decrement does not show 
any marked rise until the elastic limit has been 
greatly exceeded. Thus for crystal No. 4 (after 
the anneal, Fig. 2) a noticeable rise starts at 
about 230 g/mm? which is twice the elastic 
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limit to be expected” from the observations on 
crystal No. 6, of the same purity. For crystal 
No. 2 (Fig. 1) the decrement has only doubled 
when the stress is four times the elastic limit. 
It has not, however, been definitely shown that 
an inelastic state can be set up in a crystal of 
such low orientation'’® by the means employed 
for the higher orientation crystals. Crystals 
Nos. 3 and 5 show level decrements to very 
much higher stresses. On account of their greater 
impurity the elastic limit may well be much 
higher than for crystal No. 6. Xo statement can 
therefore be made as to how much the elastic 
limit is exceeded in these crystals. 

The observation that the decrement may re- 
main fairly constant well beyond the static 
elastic limit seems to be connected with the 
initial slow rate of transient creep previously 
mentioned. With a rapidly alternating stress it 
seems likely that sufficient plastic deformation 
to cause an appreciable hysteresis loop, and a 
corresponding increase in the decrement, does 
not have time to occur. If the presence of dis- 
locations is essential for creep and if the rate of 
creep depends on the number present,’ then 
both the rapid and slow speed experiments seem 
to point to the necessity of supposing that some 
sort of incubation period exists for generation 
of the many more dislocations needed to give 
either a rapid rise in decrement or the greatly 
increased rate of creep shown by the steep parts 
of the curves of Fig. 7. 

If, however, the rate of creep depends on the 
rate of generation of dislocations it must be 
assumed that these are generated only at a low 
rate for the first small time interval following 
an increase in stress. 

The difference found here between the two 
pairs of crystals of different purity may explain 
the behavior of the very pure crystals of refer- 
ence 1. Such pure material may have so low an 
elastic limit that it is just about reached when 


Computed on the basis of the same resolved shear 
stress. 
_ ¥F, Seitz and T. A. Read, J. App. Phys. 12, 470 (1941). 
See in particular pp. 471 and 478. 
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the stress amplitude is the lowest at which the 
measurements can be made. On this hypothesis 
the crystals would be in the inelastic state when 
mounting in the apparatus was completed, since 
strains due to handling would be unavoidable. 
The decrements observed immediately would be 
high, but would decline with time as the crystals 
annealed and approached the elastic state. The 
decrements could not be expected to stay level 
over any considerable range of stress amplitude, 
however. 

Miller’s'' experiments bear out some such sup- 
position as to the very fow elastic limit of very 
pure zinc since he obtained creep with a resolved 
shear stress equal to about one-third of that at 
the elastic limit of crystal No. 6. His creep rate 
(averaged for the first 80 hours) is roughly one- 
tenth of the steady creep rate taken from the 
last part of curve XV in Fig. 7. Such a rate would 
be well within observation with the present 
apparatus and definitely does not occur for a 
well-annealed zinc crystal of 99.99 percent purity. 
Of course if Miller’s crystal was in the inelastic 
state, due to handling, etc., both the rate of 
creep and the stress at which it occurred seem 
possible without assuming a vanishingly low 
elastic limit. The change to a much faster rate 
of creep, however, after the strain had reached 
5X10- seems to point to the onset of an in- 
elastic state at this time rather than to its being 
present earlier. 

Attempts to compare creep rates and internal 
friction have been made by Seitz and Read." If 
the mean transient creep rate taken from the 
steep parts of the curves of Fig. 7 is used, with 
a corresponding stress of 65 g/mm?, a decrement 
of about 10-* results. While this is of the right 
order of magnitude, the result must be regarded 
with some suspicion if there is any validity in 
the preceding speculations, in which it has been 
assumed that in the region of level decrement this 
rapid transient creep does not occur. 

The writers are indebted to Mr. George Kuck 
for making the spectrographic analyses. 


‘4 See Fig. 40, reference 12. 
15 Reference 12, Eqs. (4) and (5). 
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Design of Broad Band I. F. Amplifier. Part II 


RICHARD F. BauM 
Raytheon Manufacturing Company, Waltham, Massachusetts 
(Received May 13, 1946) 


This paper is a continuation of a previous article on the 
design of broad band I.F. amplifiers of the stagger-tuned 
type. The analysis is extended to unrestricted band width. 
It is shown that an exact solution is possible for either an 
oscillatory or a monotonic response, with this result: 1. The 
necessary number of stages (¢) for a given minimum attenu- 
ation in the cut-off region depends again only on the gain 
tolerance and the desired response. 2. Attentuation minima 
and maxima again appear at frequencies within the band, 
easily located and dependent only on the number of 
stages. 3. The amplifier consists of a number of pairs of 


(A) INTRODUCTION 


| iy a previous paper on the design of broad band 
1. F. ampliers,'! a method of calculating the 
circuit parameters was outlined, which leads 
to simple formulas both in the case of monotonic 
or oscillatory amplifier response. The formulas 
derived there give a very good insight into the 
influence and mutual dependence of the given 
and the calculated design quantities. This sim- 
plicity was achieved at the expense of band- 
width. It was necessary to introduce an approxi- 
mation for the expression of the detuning e. By 
replacing its exact value 


e=f/fo—fo/f 


by its approximate value 


e=2(f—fo)/fo, 


this approximation restricted the validity of all 
derived formulas to a band width not larger than 
say one-quarter the center frequency fo. 

It is obvious that a design procedure per- 
mitting the exact calculation without restriction 
of band width would be welcome. This problem 
was successfully solved by Wallman? for the case 
of a monotonic response curve. In this paper the 
calculation is extended to the oscillatory case, 
which is often preferred to the monotonic case 
for its higher gain band width product and 
steeper cut-off slope. 





1R. F. Baum, J. App. Phys. 17, 519 (1946). 
2H. Wallman, ‘Stagger: tuned I. F. amplifiers,” M.1.T. 
Report 524, February 1944. 
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circuits, which have the same figure of merit Q, and tune 
at resonant frequencies fo, disposed in geometric symmetry 
around the middle band frequency fo thus making fon 
equal to 5, fo and fo/6,, respectively. Formulas for Q, and 
6, are derived. The circuit impedances are calculated from 
a prescribed gain or from the maximum obtainable gain. 
A formula for the maximum gain band width product is 
given. A comparison is made between the performance of 
amplifiers with monotonic and oscillatory response. A 
practical example is worked out. 


Any particular desired response curve is sub- 
stantially described by 4 points. Two of them 
define:a band width within which the gain 
does not vary more than a definite prescribed 
amount. This amount is small and makes the 
gain in this region essentially constant. The two 
other points prescribe a certain minimum attenu- 
ation which has to be reached at certain fre- 
quencies outside of the pass band. The amount 
of required attenuation at these points as well 
as the admitted gain variation are given quan- 
tities for any particular design problem. In some 
cases two points may coincide with the half- 
power points; but formulas which are not re- 
stricted to the use of the half-power points are 
obviously more general and convenient. 

Wallman has shown that by arranging the 
circuits in pairs of same Q, and with resonance 
frequencies fo; and fo: symmetrically disposed 
around the middle frequency fo so as to make 
for=Sifo and fox=fo/51 a monotonic response 
curve of the form 1+(Qe)** can be obtained, 
from which the individual Q, and 6, can be 
derived by factorization. In the present paper it 
shall be shown that the most general expression 
is a polynomial in ¢ of highest power k,, = 2t, con- 
taining only even powers of e. The same line of 
reasoning is applied as in the previous paper; 
the polynomial is identified with an expression 
derived from Tschebyscheff’s polynomial, the 
roots of which have been calculated. They fur- 
nish the Q, and 6, for the oscillatory case. 

The same designations are used, as in the first 
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paper, but some of the quantities will have to 
be defined in accordance with the more rigorous 
requirements of the present calculation. Because 
of the geometric symmetry, all corresponding 
frequencies will have the middle band frequency 
as their geometric mean, for instance the upper 
and lower band limits, the two points of pre- 
scribed attenuation outside the band, pairs of 
maximum or minimum gain, etc. The geometric 
symmetry of the response curve suggests further 
the introduction of a logarithmic frequency scale. 

As the line of reasoning is very similar to that 
of the first paper, the deduction of formulas and 
the arrangement of paragraphs proceeds in the 
same sequence. In order to distinguish between 
the formulas of the first and second part the 
formulas of the second part have an additional 
index 2. 

The introduction of a reference circuit, as in 
the first part, proves of only restricted value in 
the present case as no great simplification of 
formulas can be obtained. 


(B) POSSIBLE ATTENUATION CURVES 


With the same designations as in the first part, 
the gain of a single-tuned stage is: 


Gn = 8nRn/(1+JQn€n) (2-1) 
and its absolute value 
|Gn| =gnRn/(1+Qn7€n”)). (2-2) 
The exact expression for the detuning e, is 
én =f/fon—fon/f. (2-3) 


Again a reference circuit is introduced, which 
resonates at mid-band frequency fy and for 


which* 
=f/fo—fo/f. 


A logarithmic coordinate system is to be used. 
The variables become 


8, =Ilg (f, fon); 


(2-4) 


B=lg (f/fo), (2-5) 
thus 
€, =2 sinh B,, 


¢=2 sinh £. (2-6) 


The resonance frequency fo, of the nth circuit 


* Quantities without index n refer to the reference circuit. 
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Fic. 1. Response curves of 5 stages. The center stage 
has the lowest Q and tunes at the geometric mean fre- 
quency fo. The other stages are arranged in pairs of equal 
Q, and with resonance points shifted the same distance 
above and below fo on a logarithmic frequency scale. 


shall be determined by a factor 6,: 


fon = 5nfo. (2-7) 
Its logarithm is 
AB, =lg bn, (2-8) 
which makes 
€, = 2 sinh (8—AB,). (2-9) 


In the new coordinate system the origin cor- 
responds to f/fo=1, that is, to the mid-band 
frequency. The fact that the mth circuit reso- 
nates at 6,f9 thus appears as a simple shift of its 
response curve by an amount of A@,, to the right 
\Fig. 1). The gain becomes 


(Ga| =gnRnl1+Q,74 sinh? (8—A£B,) }~', (2-10) 


which shows that the response curves are sym- 
metric with respect to the point 6=A8,. 

The index n=1 again is given to the circuit 
resonating at the low frequency end of the band 
and the index »=t to the circuit resonating at 
the upper end. The indices n=2 and n=t—1 are 
reserved for the two adjacent stages, etc. 

The gain of ¢ stages is 


a= 


\G,!| = [[ grR,[1+40Q,? sinh? (8 —Ag,,) }-}. (2-11) 


n=1 


At the mean frequency, 8 equals zero and the 
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gain becomes: 


n=t 


'Go| = I] gnRn[1+40,2 sinh? AB, }-}. 


n=1 


(2-12) 


If all stages should tune to the same frequency 
the gain would be: 
n=t 
Go| = TL gaRn- (2-13) 
n=1 
The loss in gain due to the displacement of the 
resonant frequencies is thus 


n=t 


P= T] [1+40,2 sinh? 48, J, 


n=1 


(2-14) 


and the gain can again be written in the form 


n=t 


Go => (1 l) Il £nRn. 


n=1 


(2-15) 


The part of (2-10) which is a function of fre- 
quency defines all possible response curves and 
is: 


n=t 


Dp=&"= J[[1+40Q,? sinh? (@—A8,) ], (2-16) 
1 


wed 
where d is the voltage attenuation in Nepers. 

Now the circuits may be arranged in pairs of 
the same Q, and of resonance frequencies dis- 
posed in geometric symmetry to both sides of 
the mid-band frequency.* On a logarithmic scale 
this corresponds to a shift of A, in the positive 
sense for one circuit and in the negative sense 
for the other. The response of such a pair is: 


[1+2jQ, sinh (8—A8,) ] 
<(1+27Q, sinh (8+A8,) ] 
= 1—20Q,7[.cosh (28) —cosh (2A8,,) } 
+ 740, sinh 8 cosh AB, 
=(1+4Q,? sinh? Ag, ]—4Q,? sinh? 6 
+ 74Q, sinh 8 cosh A@,, 
=A —4(Q,” sinh? 8+ 72BQ, sinh 8. 


Reintroducing ¢ from (2-6) the response of one 
pair becomes 


Dp' = [A — Q,7é |? oF 0,7 B?, 


where A and B are constants. 


(2-17) 





* This idea and the corresponding proof is contained in 
Wallman’s paper, reference 2. 
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This equation shows that the response curve 
of a pair of circuits is a polynomial of the fourth 
power of the variable e, and containing only even 
powers of e. The response of all stages is a poly- 
nomial obtained by multiplication of several such 
polynomials and can therefore possess only even 
powers of «. Its highest power ,, will be Rk», = 2t 
or twice the number of circuits. This fact is ob- 
viously not altered by the addition of one single 
circuit resonating at mid-band frequency fo, the 
response of which is quadratic in e: 


Dp=14+@e. 


(C) THE REQUIRED RESPONSE CURVE 


The result obtained checks closely with the 
one obtained in the first part of this paper, the 
only difference being that the variable e now is 
defined by its exact value (2-4) instead of its 
approximate value according to Eq. (1-6). Thus 
the same reasoning may be applied to the present 
case as to the former: 

It was shown that the required response curve 
could be represented in either of two forms: 


D,= |do| ToC) +(1+do) (2-18) 
for the oscillatory case, or 
D,* =1+2doa** (2-19) 


for the monotonic case. In these expressions do 
is the admissible gain variation within the pass 
band, and the variable a was defined by Eq. 
(1-20) with 


a=s/s,=e/eL. 


(2-20) 


This €, is that value which the variable e takes 
at the ends of the band. This definition still 
applies to the present case although e, now is 
given by 


er=f1/fo—fo/fi, (2-21) 


where f, is the upper frequency limit of the band. 
Setting 
fr=6 fo, (2-22) 
ez, becomes 
€p=6,—1/6,=2 sinh BL, (2-23) 


where 8, =lg 5, and is the equivalent of €, on 
the logarithmic frequency scale. 


The variable a thus becomes 
a=e/e;=2 sinh B/e,=sinh B/sinh By. (2-24) 
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It was shown that the points of minimum 
attenuation within the band occur at points 
uniquely determined by the number ¢ of stages 
employed and given (see Eq. (1-23)) by 


a," =cos [ (2m +1 wns 
m=Q(0, 1, 2, ---(t—1) 


=n—1 


To this corresponds a certain detuning ¢,™'", 
a certain frequency f,™'", and a certain shift 
8," on the logarithmic scale, given by 


(2-25) 


(2-26) 
(2-27) 


a,™" = e,min éL= 2 sinh B,™in €L 


B,min as lg (f= ‘fo) = lg 6,™i", 


and 


Furthermore, it was shown that it is possible 
to calculate the complex roots ao, of (2-18) and 
(2-19). Their real and imaginary parts were desig- 
nated with ar, and al,, respectively : 


” Gon = AR, + jal. (2-28) 


Their values were (Eq. (1-34) and (1-35)) for the 
oscillatory case [with added index n]}: 


ak, =cosh (c/2t) cos[(2m+1)x/2t], (2-29) 


al, =sinh (¢/2t) sin [(2m+1)x/2t], (2-30) 
where (Eq. (1-32)): 
1+d,y 
¢=cosh" : 
0 





and for the monotonic case (Eqs. (1-38) and 
(1-39)) 


aR»* = (2dy |-"?* cos [(2m+1)x/2¢], 
aly, * = [2do }-/** sin [(2m+1)x/28]. 


(2-31) 
(2-32) 


The roots corresponding to ap, are, in the 
logarithmic coordinate system, given by Bo, : 
Con = €on/€x =2 sinh Bo,,/ex. 


(2-33) 


Now as a, is complex, Bo, obviously is complex 
also: 


Bon =BRn + jBIn. 


introduced (2-33) furnishes, after 
separation of real and imaginary components: 


(2-34) 
This into 


aR, =— sinh BR, cos BI, 
€L 


(2-35) 


al, =— cosh BR, sin BI,. 
er 


(2-36) 
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These two equations can easily be solved for 
Br, and $1,. The solution is 





bh aee)==+] (=) 441], ear 
sn R,) =— ae eee ’ - 
cos ; ; r | ) 
On? On2t\2 Tr : 
cos (281) = -=4| ( ) +41] , (2-38) 
4 4 2 
where 
on? =€17>(aR,?+al,”), (2-39) 
ra? = €12(aRy? — a2). (2-40) 


By subtraction and by multiplication of (2-37) 

and (2-38) it is found that between Br, and 81, 

the following relations must exist: 
on 

—=sinh? 6r,+sin? Bip, (2-41) 


and 


+1 =cosh (28R,) cos (281,). (2-42) 


The quantity ¢, introduced here for con- 
venience obtains significance in the calculation 
of the gain. 82, is to furnish the resonance shifts 
AB, and 67, is to furnish the Q, of the individual 
circuits. 

Equations (2-35) to (2-42) apply as well to 
the monotonic as to the oscillatory case, de- 
pending upon which roots, (2-31), (2-32), or 
(2-29), (2-30), are introduced. 


(D) CALCULATION OF CIRCUIT PARAMETERS 
a. Selectivities and Resonant Frequencies 


By equalizing the roots given by Eq. (2-16) 


n=t 


Dp= TI [1+4@,? sinh? (8o,—A8,) ]=0, 


n=l 


(2-43) 


with the roots 89, as calculated in the previous 

paragraph, all circuit parameters may be cal- 

culated. By setting each factor of (2-43) equal 
to zero, its root Bon =B8R,»+j81, is found from: 
1+4(Q,2 sinh? (Bo,—A8,) =0, 

sinh (6&,+ j8In—ABn) = 7/2Qn. 

The right side of this equation is imaginary; 


thus on the left side the real component has to 
equal zero, or 


(2-44) 


AB, =BRn. (2-45) 
This leaves 
sin 61, =1/2Q,. (2-46) 
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| Equations (2-45) and (2-46) determine the 
simple relations between Q, and A@, (the selec- 
tivities and resonant frequencies) and the quan- 
tities Br, and #1,, as given by Eq. (2-37) and 
(2-38). Thus the problem is solved in principle. 
Whereas in the approximate case (treated in 
Part I) the Q, and the resonant frequencies 
depended separately on the real or imaginary 
root components @k, and alin, respectively, it is 
seen, that in the rigorous calculation they depend 
on both components in a rather involved manner. 
For this reason the result cannot be simply 
stated. Thus the importance of the approximate 
solution lies partly in the way it indicates in 
simple formulas the trend of the phenomena and 
the interconnection of the design parameters, 
whereas the final design can be based on the 
exact solution. 

For both the oscillatory and the monotonic 
response, expressions (2-39) and (2-40) can be 
simplified : 

In the oscillatory case with (2-29) and (2-30): 

€." 
0, = “ {cosh (c/t)+cos [(2m+1)x/t]}, 
(2-47) 
=e,” {sinh* (c/2t)+cos? [((2m+1)x/2¢t]}, 


r8=—{1-+e0sh (c/t) cos [(2m+1)m/t]}, (2-48) 





in the monotonic case with (2-31) and (2-32): 
on™? =er>[2do}', (2-49) 
Tn** =€77[ 2do |“! cos [(2m+1)x/t]; (2-50) 
m=n—1=0, 1, 2, ---(t—1). 


In the beginning the gain was written in terms 
of a reference circuit, which was supposed to 
resonate at mid-band frequency fo and has a 
Q, equal Q (without index). This circuit requires 
the following solutions (from (2-45) and (2-46)): 


BR, = AB, =0 
and 


sin 81, =1/2Q 
which are obtained by setting 
(2m+1)x/2t= 2/2. 
This gives, in the oscillatory case (Eqs. (2-29), 
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(2-30), (2-35), and (2-26)): 
ar=0, ar=sinh (c/2t)=1/Qer,° (2-51) 


and in the monotonic case (Eqs. (2-31), (2-32), 
(2-36), and (2-37)): 


ar*=0, 


=[2do}"'=1/Qez. (2-52) 


These relations correspond exactly to formulas 
(1-45) and (1-49) of the first part. They ob- 
viously give a very simple relationship between 
the band width e,, the admissible gain variation 
ao, and the Q of the reference circuit, which in 
the case of an uneven number of stages coincides 
with the middle stage. The corresponding value 
of « becomes in both cases 


o=o*=1/0. (2-53) 


From (2-49) furthermore (in the monotonic 
case): 


o,* =o0*=1/Q constant for all stages. 


b. The Minimum Number of Circuits 


Nothing has to be added to the corresponding 
chapter of Part I. All formulas derived there 
remain valid in the present case and this without 
restriction for the values of the quantity a,, 
appearing in formulas (1-57) to (1-64), the value 
of which is now defined by two frequencies, fini 
and fm2 (fmi1- fm2=fo?) at which the attenuation d 
should attain at least a prescribed minimum 
value d,,. Writing 


€m =fmi/fo —fo, Tents (2-54) 
ex=fr/fo—fo/fi. F 


as in (2-20), am is given by 


and 


Otm = En, /EL.- (2-55) 


c. Impedances and Gain 


The loss in gain due to the staggering of the 
resonance frequencies was (Eq. (2-14)): 


n=t 


=I [1+40,? sinh? Ag, ] 


- - T1400 I (s5 


By introduction of Eqs. (2-45), (2-46), and (2-41) 





) +sinh? a 
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this expression reduces to 
n=t 
l= [] Q,o,. (2-56) 
n=l 


The product ¢,Q, of each circuit determines its 
loss in gain (or attenuation) at mid-band fre- 
quency. This loss, of course, is the same for 
each circuit of a circuit pair for reasons of sym- 
metry, but is different for each pair. 

The product given by Eq. (2-56) thus consists 
of a series of pairs of equal factors (squares). A 
single stage resonating at mid-band frequency 
obviously does not contribute to the gain loss, 
because, according to Eq. (2-53): ¢=1/Q. 

In the case of monotonic response the product 
(2-56) is simplified because o, is the same (see 
Eq. (2-49)) for all circuit pairs: 


n=t 


I* =e," 2do}* TI] On. (2-57) 
n =] 


The expressions for the gain loss due to 
resonance staggering can be written also in 
terms of the Q of the reference circuit by intro- 
ducing o, and Q from (2-51) or (2-52). They 


become 
— Ca cos® | (2m +1)x/2t}]! 
<7 eee ek | (2-58) 
=i 


a sinh? (c /2¢) 
for the oscillatory case, and 


n=t 0, 
*= TT — 
=< 


(2-59) 


for the monotonic case. 

The last equation shows that in the monotonic 
case the gain loss /* increases for each stage in 
proportion to the increase of Q, over the Q of 
the reference circuit. But according to Eq. (2-58) 
in the oscillatory case this loss increases more 
rapidly and differs from circuit pair to circuit 
pair. 

The loss in gain due to the staggering of 
resonance frequencies is, of course, totally or 
partly compensated by the increase in stage gain 
due to the increase in Q, as will be apparent from 
the following gain formulas. 

The total gain according to Eq. (2-15) was 


n=t 


Go| = (1/1) II 2.R,, 


n=1 
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and with (2-56) 


n=t 
Go| = [] gnRn/Qnon (2-60) 
n=l 
n=t 
ws Il Rn/ CrwonFn- (2-61) 
n=l 


As any circuit pair resonates at frequencies the 
product of which equals the square of the mid- 
band frequency, one may replace, in the product, 
won by wo. Furthermore, considering only the 
case of stages with the same capacities C,=C, 
and the same transconductance g,=g, the gain 


becomes: 
g t n=t 1 
Go -| =| II —. (2-62) 


Wo n=1 On 


Writing this expression for just one pair of 
circuits it is easily seen that in the monotonic 
case the average gain per circuit pair equals the 
gain of the reference (or middle) circuit. 

The average gain per stage is 


g n=t 1 l/t 
El = A II — | : (2-63) 
WOOL n=0 Tn 


The total band width (with Eqs. (2-22) and 
(2-23)) is 
’ fo 
BW = bi fo— = fer, (2-64) 


L 


thus the average gain band with product becomes 


ro g ez{ == 1 1/t 
Go pw-(+.)-<| II ~| . (2-65) 
rC 2 n=1 On 


For the monotonic case the root equals Q and 
with (2-52) the expression becomes 


|Go*| BW = (=) -3[2dy }'/** = (2-66) 
rC 

which is exactly the same as the expression 
derived for the approximate calculation (see 

Eq. (1-78)). 
The factor in parenthesis again is recognized 
as the theoretical maximum gain band width 
product and again an efficiency factor may be 


defined by: 
€L n=t 1 1/t 
| II —| ; (2-67) 


n=1 Tn 
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Fic. 2. Efficiency of amplifiers with oscillatory and 
monotomic response. 


or with (2-47): 


n=t 
u=2T] {sinh? (c/22) 


n=1 


(E) COMPARISON OF MONOTONIC AND 
OSCILLATORY RESPONSE 


The deduced formulas permit the comparison 
of the performance of amplifiers of monotonic 
and oscillatory response. 

For a given response characteristic, that is for 
a given maximum gain tolerance do within the 
pass band and a given minimum attenuation d,, 
above a given frequency, a monotonic amplifier 
requires more stages than the oscillatory one. 
The ratio of their number is, with Eqs. (1-58) 
and (1-63): 


t* h* cosh" am 


t ke 





(2-70) 
lg am 


The first factor almost equals one, so that ap- 


+cos? [(2m+1)x/2t]}—"/* (2-68) : 
proximately 

and in the monotonic case: 

t* cosh an 
u* =3[ do}, (2-69) — = —__—__, (2-71) 

t lg am 

the latter approaching rapidly } with increasing il ; 

the number of stages. thus if the minimum attenuation is defined 


The importance of the quantity ¢, in all gain 
calculations is apparent. With reference to Eq. 
(2-62) one can say that o, defines the relative 
contribution of each individual stage to the 
total gain. In order to facilitate the gain cal- 
culation, the value of u and u* is plotted in Fig. 2 
as a function of ¢ for dp =.05, d)=.1, and dyp=.15. 


anywhere between a@mn=Qm*=1.1 and 2, this 
ratio varies between 4.5 and 1.9. 

For a given gain tolerance within the pass 
band, and with the same number of stages, the 
monotonic amplifier has a much flatter cut-off 
characteristic than the oscillatory one. Thus by 
setting ¢/t*=1 in Eq. (2-71) the frequencies at 


TABLE I: Calculation of the Q, and resonance frequencies. 


























stage m (2m+1)x/t=0 cos 6 cosh ¢c/t+cos @ 
land5 | 0 and 4 36° and 324° .809 2.006 
2 and 4 1 and 3 108° and 252° — .309 .888 
3 2 180° ~1 
Eq. 
stage cosh c/2t-cos 6 on on? /4 (o,2/4)? tr2/2 
land 5 .97 es Ads .0305 .343 
2 and 4 —.37 .310 .0776 .006 .109 
3 
Eq. (2-47) (2-48) 
stage o cosh 28R, cos 2BT» 2BRn 2BT» AB, 
land5 | .837 1.345 .995 .197 3.7" .398 
2and4 | 557 1.136 980 516 11.3° 258 
3 262 0 
Eq. (2-53) (2-37) (2-38) (2-45) 
stage sin BT), On bn fon 
land5 | .050 10 1.49 17.9 and 8.05 mc 
2 and 4 .0984 5.07 1.29 15.5 and 9.3 mc 
3 3.82 12 mc 
Eq. (2-46) (2-8) (2-7) 
(2-53) 
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which a given attenuation is reached are given by 


Cm * = m+ (Ay? —1)?. (2-72) 
Within a range of 1.1<a,,<2 the corresponding 
am* lies between 1.55 and 3.7. 

For a given minimum attenuation and a given 


number of stages the gain variation of the 


monotonic amplifier is obviously larger. 

The average gain band width product, with a 
given number of stages and a given gain vari- 
ation, is larger for the oscillatory amplifier by a 
factor: (Eqs. (2-68) and (2-69)) 


Go -BW be 
> =—, (2-73) 
n* -BW u* 
The value of » and yu* was plotted in Fig. 2. 
Whereas the value of u* is somewhat smaller 
than 3 the value of yu is smaller than 1. The 
limits of uw can be seen to be 


(2m+1)r 


2t 


u—} sinh (c/2t) for sinh (c/2t)>>cos 


and 


p—2-''* for sinh (c/2t)-0. 

The first case corresponds to a very small 
gain tolerance and a small number of stages. 
It falls to values below one-half, but is still larger 
than the corresponding u* of the monotonic 
amplifier. 


The other case corresponds to a large gain 











“3 -4 =- “2-2 604482 3 4 


Anal 


Fic. 3. The resultant response curve shows an oscillatory 
behavior within the pass band 287, defined by a given gain 
tolerance dy. The attenuation exceeds a prescribed mini- 
mum value d,, outside a given range 28». 


average efficiency increases with an increasing 
number of stages; the total gain band width 
product tends towards 3, whereas in the mono- 
tonic case the average gain band width product 
tends toward 3. 

From the standpoint of amplitude response, 
which was considered exclusively, the oscillatory 
amplifier is superior to the monotonic one in 


tolerance and a large number of stages. The every respect. 
TABLE II 
(2m+1)9r @ @ 
m ” 2 ve 2 €,™u By ) 
O and 4 18° and 162° 951 .794 387 1.473 8.16 and 17.7 mc 
land 3 54° and 126° 588 .490 .242 1.274 9.42 and 15.3 mc 
2 90° 0 0 0 | 12 me 
formula (2-25) (2-26) (2-26) (2-27) 
TABLE III. 
mr =) 
” t cos/ y €n 4 B»,™x bm" * ‘x 
land 4 36° and 144 .809 .675 331 1.392 8.64 and 16.7 mc 
2 and 3 72° and 108° 309 25 .129 1.138 10.55 and 13.7 me 
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(F) EXAMPLE 


A stagger tuned |.F. amplifier is to be designed 
having a band width from 8 to 18 mc, and having 
the maximum gain obtainable with 6AKS5 tubes, 
with a gain variation of .1 Nepers within the 
pass band, with a power attenuation of at least 
30 db at 22 me and the corresponding lower fre- 
quency. An oscillatory response is desired. The 
given quantities are: 

dy=.1 Neper, the admissible gain variation. 
é~ = 1000, the required power attenuation. 
f;=18 mc and 8 me, respectively, the upper 
and lower band limits. 
fo =(8.18)?=12 me,: the geometric mean fre- 
quency, 


12° 


fr. =22 —=6.55 mc 
22 


the frequencies at which the power drops to the 
required level. To f;, and f,, correspond detunings: 


é,= .83 
and 


€m = 1.284 
which correspond on the logarithmic scale to 


BL = .407 
and 
Bm = .607 
(see Fig. 3). Thus 
cosh~ a,, = 1. 


Qm = €m /€,= 1.542 and 


The minimum number of stages (from (1-58) 
and (1-59)): 


2t=h/cosh™ a, 
cosh h = [é4" —(1+do) ]/do 
=10000 


h=9.8. 
=9 8 


In order to obtain the prescribed response curve, 


at least 5 stages are required. Calculation of the 


hyperbolic expressions: 


1 +d, 
c=cosh™ - 7 —=cosh™ 11 =3.09, 


0 
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. 3.09 
C/ 2t ao = .309, 
10 


c/t=.618, 
cosh (c/2t) = 1.048, 
sinh (c/2t)= .314, 


cosh (c/t) =1.197. 


The resonance frequencies fo, and the Q, of the 
circuits have been calculated in Table |. The 
result has been plotted on a logarithmic scale in 
Fig. 1. By addition of the 5 individual response 
curves the total amplifier response of Fig. 3 is 
obtained which shows the required oscillatory 
behavior. The location of minimum attenuation 
is shown in Table II. 

The points of maximum attenuation can be 
calculated in an analogous way from: 


mr 
Q,,'"** = cos (“*). m=1,2,3---(t-—1 ), (2-74) 
t 


which follows from inspection of Fig. 3 of Part I, 
(Table III). 

The corresponding maximum and minimum 
values of attenuation within the band turn out 
to be 5.46 and 5.54. The former checks with the 
loss in gain if calculated from Eq. (2-58). The 
difference of .08 checks with the admissible gain 
variation, as logo (1.1)?=.08. (All curves are 
power response curves.) The minimum attenu- 
ation of 30 db=3 is reached within the pre- 
scribed frequency limit 8,,=+.607 with an 
attenuation of 8.57. 

Assuming a transconductance of g = 5000 mhos 
and a minimum capacity of 12uyf for each stage, 
the theoretical maximum band width 
product is: 


gain 


9 /rC=133 me. 


The efficiency factor is. calculated from Eq. 
(2-66) with 
833 _ a 
a“ t (.837? XK .557* K .262)? =.73. 


The average gain band width product is thus: 
Go| BW =133X.73 =98.5 me. 


With a band width of 10 me the average gain 
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becomes 
6, = 98.5 ‘10=9.85, 


and the total gain: 
|Go| =9.85°= 92,400. 


This compares very favorably with an amplifier 
of monotonic response whose best efficiency 
factor is .5. This gives a maximum gain for 5 
stages (with 10 mc band width) of only 13,000. 


CONCLUSION 


It has been shown that the stagger tuned am- 
plifier problem with oscillatory response can be 
exactly solved mathematically. The solution is 
obtained by the same means employed for the 
solution of the approximate case (with restricted 
band width), that is—by calculation of the roots 
of either a parabolic or a Tschebyscheff poly- 
nomial. 





The Numerical Solution of Laplace’s Equation in Composite Rectangular Areas 


Max M. Frocut* 
Illinois Institute of Technology, Chicago, Illinois 


(Received January 29, 1946) 


An analytical method is developed which yields key 
values of a harmonic function U in rectangular figures in 
terms of certain simple combinations of known boundary 
values, which we call Laplacian perimeters and denote by 
P. From the key values all the remaining interior values 
can be directly and rapidly calculated. The basic ideas 
consist of the introduction of four elementary figures of 
variable lengths: 3(2+6), 3(3+6), 4(2+46), and 4(3+5), 
where 0£ 6< 1, for which the key values are found alge- 
braically, Eqs. (9.1) to (9.8) inclusive, and of combining 
analytically these elementary cases to obtain continuous 
solutions for rectangles (3Xm’) and (4Xn’) for any 
number n’ greater than 2. The general expression for any 
key value (ui )m xn in a rectangle (mXn’) is 


1. INTRODUCTION 


UMERICAL or approximate solutions of 
Laplace’s equation 


#U/dx?+aU/day? =0 (1.1) 


are generally based on the fact that the value of 
’ the function U at an arbitrary point O of the 
region R to which the function U applies, Fig. 1, 
must satisfy approximately the general equationt 


1 1 Ua Us ) 
a 
ac bd a(a+c) b(b+d) 


U- Ua (1.2) 
+—_—_—_-+——_- (a) 
c(c+a) d(d+b) 
or 
Uy =RattatRytn+ku-t+kata, (b) 
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(Ui)myen? = CiP1 + CoP2+ +++ CePi +--+ Crepe, 

in which C;, C2---C, are constant coefficients. Analytical 
relations are established between the coefficients for the 
last key value u, in a rectangle (mXmn), where n is an 
integer, and the coefficients for the next key value ue4: 
in the rectangle m(n+6) or m(n+6+1). Specific values 
of the coefficients are presented in tabular form for rec- 
tangles (3Xn’) and (4Xn’) with the aid of which all key 
values in such rectangles can be rapidly calculated. The 
tables also permit the calculation of any one interior key 
value by itself. A numerical example and further approxi- 
mate formulas are given, and applications to composite 
rectangular areas such as angles, channels, etc., are dis- 
cussed. 


in which ug, us, u-, and ug are the values of the 
harmonic function U at four neighboring points, 
A, B, C, and D, respectively, with coordinates 
as shown in Fig. 1. Equation (1.2) will be referred 
to as the four-point influence equation, or more 
briefly the four-point equation. 

The procedure in numerical solutions generally 
consists of drawing a rectangular, equally spaced 
network distant h apart, and adjusting the 
values at each point of intersection of the net- 
work lines (network points) until they satisfy the 
four-point equation. When that is attained, an 
approximate or numerical solution is said to have 
been found. The distance h between the network 


* Professor of Mechanics. 
+ For a derivation see references 1 and 3 in bibliography 
at end of paper. 
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x 
A— 
(a, 0) 


D (0,-d) 
Fic. 1. 


lines is called the lattice unit. The points A, B, 
€, and D may be located on the network lines, 
or on lines inclined 45° to them. In the first case 
we speak of normal neighbors, in the second case 
of diagonal neighbors. In the treatment which 
follows, both types are used. 

Most of the interior network points will have 
equally spaced neighbors. For such points (1.2) 
reduces to 


4u, =Ugtuptuetug. (1.3) 


Equation (1.3) is generally called the Liebmann 
formula. 

Another case of special interest is the case in 
which a=), and c=d, Fig. 1. For this case (1.2) 
reduces to 





C(Uatuy) a(u-t+tUa) 
= 4. . 
2(a+c) 2(a+c) 


Uo 


(1.4) 


Assuming a <c and a/c =6, this may be written as 





(Uatuy) 5(u.+ua) 
Up = aa , (1.5) 
2(6+1) 2(6+1) 
where 6 is a fraction less than unity, or 
Uo = Ky(uatus)+Ko(ue+uz), (1.6) 


where K, and Ke are constant coefficients given 
bv 


K,=1/2(6+1), (a) (1.7) 
K,=6/2(6+1). (b) |<” 
When 6=c=d, a<b, and a/b=6, we have 
Ua bu. 5(up+Ua) 
Uo 1.8) 





= + + 
(14+6)? (1+6)? 2(1+8) 
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The process of adjustment which makes the 
values of a harmonic function agree with the 
four-point influence equation at each network 
point may be called a process of harmonization. 
The most commonly used method for harmoniza- 
tion is that of iteration’ in which, arbitrarily 
assumed or judiciously selected, initial values 
are successively improved until practical con- 
vergence is obtained. Another procedure which 
has proved effective is that of relaxation.? More 
recently harmonization in rectangular areas has 
been effected by a solution of a special system of 
difference equations.?* 

The object of this paper is to develop an 
analytical process of harmonization, which yields 
key values of the harmonic function U in rec- 
tangular figures in terms of certain simple com- 
binations of the given boundary values. From 
the key values all, except two or three, of the 
remaining interior values can be calculated di- 
rectly by the simple Liebmann formula (1.3). 
For two or three points (1.8) may be required. 
The present treatment is believed to result in 
simple derivations, a relatively large scope, and 
a material reduction in the labor required to 
effect a solution. The method is directly appli- 
cable to composite rectangular figures, such as 
angles, channels, T-figures, Z-figures, etc. The 
proposed process has been used in numerical 
cases and found to be rapid and accurate. 

The rectangles treated have an arbitrary width 
m and length n’. The width is always represented 
by an integer. Specific solutions are obtained for 
rectangles (3Xn’) and (4Xn’) for all positive 
values of n’ between 2 and 14. The results are 
summarized in tables at the end of the paper. 
The tables are applicable to any value of n 
greater than 14. 


Key Values. Definition 


A set of interior values of U, so located in the 
region R to which the function U applies that, 
from these values, the values of U at all other 
network points can be calculated by the Lieb- 
mann or four-point formulas, applied either 
normally or diagonally, is defined as a set of key 
values. In a square (3X3), Fig. 2(b), any one 
of the four interior values 4, uo, Us, and uy is a 
key value for the whole region. Thus if u; be 
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found from the diagonal 
neighbors, and us, v4 from the normal neighbors. 

Similarly in a rectangle (4X 3), Fig. 3(b), either 
u, or u,is a key value for the whole region. Thus 
if uw, be known, uw; and us can be found from 
diagonal neighbors, and wo, uy, us from normal 
neighbors. 

A similar condition exists in a (4X4) square, 
Fig. 5(b), in which m1, the value at the center of 
the square, defines the remaining eight interior 


known, us can be 


values. In all three cases one key value is suf- 
ficient to obtain a complete numerical solution of 
Laplace's equation. 

In rectangles (mXn), m=3, 4, n an integer, 
the number of key points & is given by 


n-odd, 


, n-even. 


k=(n—2)/2, 
k=(n—1)/2 


In rectangles (m Xn’), n’=n+6, Of 6< 1, 


k=n 2, n-even, 
k=(n—1) 2, n-odd. 


The key points are numbered 1, 2, 3---& and 
the corresponding key values are denoted by 4, 
Uo, Ug** Ux, respectively. 

An interior, or unknown, value will always be 
denoted by u, with a suitable subscript relating 
it to a given point. Known boundary values will 
be denoted by a, 6, c---, ete. 


2. KEY VALUES FOR INCOMPLETE AND FULL 
(3X3) RECTANGLES, FIG. 2 


(a) A rectangle 3(2+ 6), in which 6 is a fraction 
less than unity, will be spoken of as an incom- 
plete (3X3) square, Fig. 2(a). 

In order to find the value u; we use (1.6) and 
(1.3). Thus, 


u; = K,(a,+a2)+K2(a3+4s), (a) 


















































(2.1) 
43 = Uy +0; +62+C3. (b) 
a, ™ - bn as a, 
Qa 
6 s C M 1% 
4%, a . 
\ 7 
- - a 
4| “4 b 4% |%, ; 
G St G, 
Sg 7 G 
(a) (b) 
Fic. 2. 
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Upon solution, eliminating K, and Kg by (1.7), 
we have 


(84-76)uy = 4(a,+a2) + 6(4a3+¢,+6¢2+¢3). (2.2) 


Putting 6 equal unity we obtain 


3 3 
15u,=4 > a;t+>d c;. 


i=1 i=] 


(2.3) 


Equation (2.3) gives the key value for a full 
(33) square. 


3. KEY VALUES FOR INCOMPLETE AND FULL 
(3X4) RECTANGLES, FIG. 3 


(a) A rectangle 3(3+6), Fig. 3(a), will be 
spoken of as an incomplete (3X4) rectangle. 
Equations (1.3) and (1.6) give the key value 4, 
for this case. Thus 


4u;=u3+u5,+d,+de, (a) | 
4u3=u,+a,+0,;4+c, (b) +(3.1) 
us=Kj(a2+62)+K2(uitei). (c)} 


Solving, and eliminating K, and Ke, we obtain 


(“== 


Jur = (artbrbertads +44) 
1+6 


2 
+(—) ar+b.+ 40) (3.2) 
1+6 


Putting 6 equal unity, we obtain the key value 
for a full rectangle (34), Fig. 3(b). Thus 


144, =a,+42+2¢,+4(d;+d2)+6,+)2. (3.3) 


By putting 6 equal zero, changing the notation 
of the boundary values to agree with that of 
Fig. 2(b), and observing that when 6 becomes 
zero d2=b2=a,; we obtain (2.3), ie., the ke 
value for a full (33) square. 
























































a, 44> - oh a a, 
h 
b 4% t y) My 1% |% 4 
/ 6, 
’ Ms 
“, | 4s a, Us; us 
a a. 
/ C a C 2 
(a) (b) 
Fic. 3. 
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4. KEY VALUES FOR INCOMPLETE AND FULL 
(4X3) RECTANGLES, FIG. 4 


(a) A rectangle 4(2+6), Fig. 4(a), will be 
spoken of as an incomplete (4X3) rectangle. 
























































4 md a "6 5A a, 4, 
4, u; u, d 
a1 \. d 
C SE 4, & 
[ 4, J d, 
4 “ ly 
S\N : 
a a 
? 
% a, 4, 
(a) (6) 
Fic. + 


To find the key value 4; we again use (1.6) and 
(1.3) and obtain the following three equations: 


uy; = K,(d,+d2)+Ko(u3t+us), (a) ' 
4u3=u,+a,+6,+¢, (b) (4.1) 
4u5 =u, +de+b.+c). (Cc): 


Upon solution we have 


4+ 36 2 
( us =“(dtds) 
6 6 





+3(aita2+bi+b2)+e,. (4.2) 
(b) Putting 6 equal to unity in (4.2) we obtain 
7) = 3 (di: +d2+b1+b2) +2(di+d2) +c). (4.3) 


This is the key value for a full (4X3) rectangle, 
Fig. 4(b), which agrees with (3.3) giving the key 
value for a full (3X4) rectangle at the same 
point. 


5. KEY VALUES FOR INCOMPLETE AND FULL 
(4X4) SQUARE, FIG. 5 


(a) A rectangle 4(3+4) will be spoken of as an 
incomplete (4X4) square, Fig. 5(a). There are 
three ways to obtain an expression for the key 
value “;. The simplest procedure is to express “4; 
in terms of us, 43, Us, and us and then to express 
these quantities in terms of ;. Using (1.3) and 
(1.6) we obtain the following five expressions: 
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du; tuetugtugtus, (a) 
4u3;= Uy +de+CotCz, (b) 
4u,=uUyt+a3+c3+C4, (c) 7(5.1) 
U2 = Ky(a,+e1') + Ko(ui+ce), (d) 
us = K,(ay+c,"") + Ko(ui+c,). (e) 


Upon solution we have 





7+56 
( Jur Hastastertdertey 
1+6 





1 
~~ (aitast+er' +1") 
(1+6) 





(Co+c4). (5.2) 
(1+6) 


(b) Putting 6 equal unity and observing that 
then c;'=c,’’=c, we obtain 


4 4 
6u,=4 > aitd ci. (5.3) 
ixl i=1 


This gives the key value for a full (44) 
square. 

(c) Equation (5.3) also yields the key value 
for a rectangle (4X3). Putting 6 equal zero, 
changing the notation of the boundary values to 






























































4 
4 pA "lee nt ro a: & a, 
a 
U4, . U4, |, |4, 
7) ¢ 4g |\% 1% 
G ~ C a . C 
%, ee 
¢ uu, |4e | us 
&, 
‘i < 9 7*| ¥4, 
J a 
4 4 & 
(a) (6) 
Fic. 5. 


agree with that of Fig. 4(b), and observing that 
when 6=0, c;’=a,;=d,, and c,;"’=a4=d2, we 
obtain (4.3) for a full (4X3) rectangle. 


6. LAPLACIAN PERIMETERS 


In the further development of the method 
under consideration, extensive use is made of 
the key values for the basic figures treated in the 
preceding sections. By combining the results for 
the elementary rectangles, solutions are built up 
for longer rectangles. It is therefore convenient 
to simplify the forms of the expressions of the 
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key values already derived. Specifically we shall 
write each key value in the form 


au=P+u’, (a) 
or, (6.1) 


au=P+u'+u", (b) j 


in which @ is a constant, P is a quantity deter- 
mined by the known boundary values, and wu’, u, 
u’’ are consecutive key values. The quantity P 
in (6.1) will, for lack of a better term, be referred 
to as a Laplacian perimeter or more briefly 
perimeter. Since we have eight basic rectangles, 
we shall need eight Laplacian perimeters: four 
for the incomplete rectangles and four for the 
full figures. We shall denote these perimeters by 
the generic symbols p and P, for the incomplete 
and full rectangles, respectively. As already 
stated, these perimeters are to contain only 
known boundary values. It will be seen later 
that, in each configuration, one or at the most 
two, boundary values will be characteristically 
missing in comparison with the full expression 
for the key value. Thus, incomplete and full 
(33) squares will always be used in such a 
manner that c; will be unknown, and in the 
full (34) rectangle both }; and. b2 or at least 
b,, will be unknown. In the definitions of the 
eight Laplacian perimeters that follow, we are 
guided entirely by such considerations. 


7. DEFINITIONS OF LAPLACIAN PERIMETERS 


Notation : The notation Pn or (P) men denotes 
a Laplacian perimeter in a rectangle (m Xn). 


Full rectangles—P Incomplete rectangles—p 


Psys—P, P3x3x—P , 
P3.4—P", P3xs— pp”, 
Pyxs—P’, Pixs—?P’, 
Poa—P, Paxa—P. 


Definitions : 
1. Incomplete (3X3) square, Fig. 2(a), 


4 
p =—(a, +42) +4a3+¢2+€3. (7.1) 
6 
2. Full (3X3) square, Fig. 2(b), 
7 3 
P=4 ¥ ajytceo+c3. (7.2) 
i=] 


3. Incomplete (3X4) rectangle, Fig. 3(a), 
bp” =a, +¢,+4(d,+d2) 


2 
+(—=)tartbet ier) (7.3) 
1+6 


4. Full (3X4) rectangle, Fig. 3(b), 
P” = (a,+a2+2c,)+4(d;+d2)+known b's. (7.4) 


5. Incomplete (4X3) rectangle, Fig. 4(a), 
2 
p’= (ai tas thi tbs) + (di tds). (7.5) 


6. Full (4X3) rectangle, Fig. 4(b), 





P’=}(ait+a2+b:+62)+2(ditds). (7.6) 
7. Incomplete (4X4) square, Fig. 5(a), 
p= }(de+a3+cC2+c4) 
1 5(Co+c,4) 
+(—)etarte’te")+— ~ Ba 
1+6 1+6 
8. Full (44) square, Fig. 5(b), 
4 
P=} > a;+> known c’s. (7.8) 


i=l 


8. RELATION BETWEEN P,,,.,, AND Pry n’ 


(a) In all figures treated in this paper there 
are as many Laplacian perimeters as there are 
key points, or key values. 

(b) The numeral subscript in the perimeter 
refers to the key point of the same number. 
Thus P; refers to key point (4) of a key value uj. 

(c) The notation (Pi) mxn denotes a perimeter 
at a point (z) in a rectangle (m Xn). 

(d) When the length of a rectangular figure is 
increased by \, only the last perimeter P; in the 
rectangle (4X) changes. In rectangles (3 Xn) 
and (4Xmn) the relation between the old and 


a & 





Fic. 6. 
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Uy, 
(a) d, 2 
leA-sd —sl le-5A 
az & z 4, ° 
a qd, 
G Cell 42 
N q, 
4 4 %* 4 4 
(4) (c) 
Fic. 7. 
new perimeters is 
(Px)mxn = (Px) mxn? +Ue+15 (8.1) 


k referring to the figure (m Xn). 

Consider for concreteness the rectangle (4 X 10), 
Fig. 6. In the last square having its center at 
point (4), the Laplacian perimeter (P4) 4x10 is by 
definition, (7.8), 

4 


(Ps) 4x10 = 3 )e Qitcitcety, 
i=1 
since ¢;=43 is unknown. If the length »=10 is 
increased by \ 


M. 


t 


a:+CotCa, 


— 
(Ps) sc104a) = 


since now both c,;=us5, as well as c3=43, are 
unknown, so that 


(Ps) 4x10 = (Ps) scoany + Us. 


We note, however, that P:, P2, and P; remain 
unchanged, or in general 


(Pi) mxn = (Pi) mxn’s 


t<k. 


for 


(8.2) 


4 eh 


di, 





(a) 
Fic, 8. 
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9. BASIC KEY VALUES IN TERMS OF P AND p 
1. Incomplete (3X3) square, Fig. 2(a). Eq. 
(2.2) becomes 
&i=pteci, (a) 
& = (8+78)/6. (b) 


2. Full (3X3) square, Fig. 2(b). Eq. (2.3) 


becomes 


(9.1) 


15u,=P+c,. (9.2) 


3. Incomplete (3X4) rectangle, Fig. 3(a). Eq. 
(3.2) becomes 


au, =p" +b, (a) 
a” =(15+135)/(1+6).  (b) 


4. Full (3X4) rectangle, Fig. 3(b). Eq. (3.3) 
becomes 


(9.3) 


14u,=P”’ +> unknown b's. (9.4) 


5. Incomplete (4X3) rectangle, Fig. 4(a). Eq. 
(4.2) becomes 

a’u=p'+c1, (a) 

a’ = (4436) /6. (b) 

6. Full (4X3) rectangle, Fig. 4(b). Eq. (4.3) 


becomes 


(9.5) 


7u,=P’+¢,. (9.6) 


7. Incomplete (4X4) square, Fig. 
(5.2) becomes 


5(a). Eq. 


au, =pt+cs, (a) 


a=(7+58)/(1+8). (bf? 

8. Full (4X4) square, Fig. 5(b). Eq. (5.3) 
becomes 

6u,=P+> unknown c’s. (9.8) 


10. GENERAL METHOD 


By combining the solutions for the values of 
U for the basic cases it is now possible to obtain 
continuous solutions for the key values in rec- 
tangles (m Xn’), where n’ is any positive number 
greater than 2 and m=3, 4. 

The general method rests on the following 
simple proposition : If a solution for the last key 
value u, is available for a rectangle (m Xn), n 
any positive integer, it is possible to extend this 
solution to wus, in rectangles m(n+6), or 
m(n+1+6). The range of 6 is O<6<1. As a 
matter of convenience it is best to start with n 
equal to an even positive integer. 


~I 
w 
wn 








Any rectangle (m Xn’) may be viewed as a 
rectangle (m Xn), n the nearest even integer, 
overlapped by an incomplete (#3) or (mm xX 4). 
Thus, a rectangle 4(6+6) may be looked upon 
as a (4X6) overlapped by an incomplete (4X3), 
and a rectangle 4(7 +46) may similarly be viewed 
as a (4X6) overlapped by an incomplete (4X 4), 
Figs. 7 and 8. Now, in the rectangle (m Xn) the 
last key value u, may be written as 


uy, = C,P\+CoP.+---C.P, (10.1) 


in which the C’s are constant coefficients and the 
P’s are perimeters for the rectangle (m Xn). 
Remembering the relation between P,,., and 
Prnsent, (8.1, 8.2) we have 


uu; = CiP,:+C.oP2+ os Cy. Pi t+ Cuteness, (10.2) 


in which the C’s are the same as (10.1) but the 
P’s now refer to the rectangle (m Xn’). Also 
from the incomplete rectangles (m X 3) or (mX 4), 
Section 9, we have 


OU je = Presi tur. (10.3) 
Solving (10.2) and (10.3) we obtain 
Mpa = (a— Cy)" (C Pi + CoP 24+ +++ CPx) 
+ (a— Cy) Dest. (10.4) 


We thus have a general relation between the 
coefficients for the last key value u, in the 
rectangle (m Xn) and the coefficients for the 
next key value “,,; in the rectangle m(n+5), or 
m(n+1-+-5). Specifically 


Crai=(a—C,), (a) } 


and for all other coefficients (10.5) 


(Ca)mxn’ = (Ci) mxn(a—Cy)—. (b)} 


In using these formulas, it must be remembered 
that there are four a’s, which depend on whether 
m=3 or m=4, and on whether n’ differs from 
n by less than unity or by more than unity. 
Thus, for rectangles 4(n+1+6) and 4(n+6) we 
use a and a’, respectively, which are given by 
(9.7b) and (9.5b). Similarly for rectangles 
3(n+1+6) and 3(m+6) we use a” and &, respec- 
tively, given by (9.3b) and (9.1b). 

The determination of u,.; reduces the rec- 
tangle (m Xn’) to that of a rectangle (m Xn) for 
which a solution for all key values is assumed 
to be known. 
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Hence, starting with the expressions for the 
key value “; in a (4X4) square, (9.8), one can 
successively write down the formulas for the 
last key values in all subsequent rectangles 
(4Xn’). Similarly, starting with the expression 
for mu, in the rectangle (34), (9.4), we can 
write the formulas for all rectangles (3 Xn’). 


11. SUCCESSIVE FACTORS (a—C;,)7'! 


For brevity we simplify our notation and write 


(a—C,.);"=F,, (11.1) 
and 
CF, C, a) mines) = (F, C, x) ars. (11.2) 
We then have 
Fy45= (C2) 443 =Laess—(C1)s | i (11.3) 
F5i8= (C2)548=[Lasis—(Ci)s ], (11.4) 


Fois=(Cs)648= Laois — (C2)6 }- 
=(ao.;—F,)—', (11.5) 
Fy.5= (C3) 748 = Lasis—(Co2)¢6 
=(a3is—Fs)“!, (11.6) 
Fis=(C3)sis = Laess—(C3)s |“ 
=(ae.3—Fs)“. (11.7) 
In general, for m an even integer 
Fiis=(Qex8—F,)“, (a) 
Fisiss=(azis—F,). =(b) 


(11.8) 


For example, to find F45,3s) we write 


F 454.8) = (@¢4¢34.8) — Fays) 


= [as¢a48) —(Crdaxa]. = (11.3) 


Substituting 
434.8) = (7 +56) /(1+5), (9.7b) 
and 
(Ci)sxs=6) (9.8) 
we have 
F 4548) = (41+296)/(1+6)6. (11.9) 


We can now write down the expression for 


(2) 4¢548). Thus by (10.4) 
(to) 4¢548) = F4¢548) (Ci) aaP it Facsisyp2o. (11.10) 


Calculating several values of F it is found 
that these factors converge rapidly, and that 
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TABLE I. 


C3 Cy Cs Ce Cr 
446 (Ci) 4 Fass F445 
5+6 (C1)4F 5.8 F3.5 
6+6 (Ci)sFeFo+s FeF 6.5 Fes5 
7+6 (Ci) a Fo F745 FeFr35 Fr48 
8+6. (Ci)sFeFsFsis Fe FsFs.5 Fs Fs F3.5 
9+6 (Ci)4F Fs Fo.5 Fe FsFo45 Fs F438 Fyo..5 
10+65 (Ci)sFoFsFioF i045 FeFsFioF 1045 Fs FioF 1045 FyioF 1045 Fy048 
114+6 (Ci)sFeFsFioF iss FoF sFioF i438 FsFyoF i145 FyoF iiss F345 


12+6 (C1)sF 6 FsFioF 12F 12.5 


FoF sFyoF 12F 1245 
13+65 (Ci)sFeFsFioFi2Fiss8 


FoF sFioF12F 1348 


FsFioF 12F 12.5 
FsFyoF 12F 1345 


FioF 12F 1245 


Fy2F 1245 Fio4s 
FyoF 12F 1348 


Fy2F 1345 F345 


14+6 (Ci)sFoFsPioFi2FisFisss FeFsFioFi2FisFiass = FsFioFi2FisFisss) | FioPioFisFisss  FioFisPises  FisFisas Pisses 


15+6 (CidsFeFsFioFi2FisFisss FeFsFioFi2FisFisss FsFioF12F 14F 1548 


they become essentially constant at n=6 when 

m=4, and at n=4 when m=3. Thus, to five 
decimal places, 

Faenss) = Fagas), modd; (a) | 

Fgenas) = Fass), neven; (b) | 

\ 

" f 

F148) = Faisis), nodd; (c)| 


+ | 
F3¢na8) = Faieasy, meven. (d) 


It is easy to verify (11.11) for integral values of 
n, i.e., 6=0. For example 


F 4, oe Fs, 7= .0000031 ’ 


Fy 1.6 Fyx-3 = .0000042 9 
F 3. yaaa F3,. 5=> .000001 65 ’ 
Fs. a F 3, 6é= .000001 89. 


But by (11.8) for 2 an even or odd integer 
Fysoss— Pros =(a— Fyse)—(a— Fy). 


Since the a’s are the same, and for integral values 
of n, Fri2=F,, for n>6, m=4, and n>4, m=3 
it follows that the relations given in (11.11) are 
correct to five decimal places. 


12. TABLE OF MASTER COEFFICIENTS FOR THE 
LAST KEY VALUE u;,,; 


From (10.5) it is possible to write down all 
coefficients from C,; to C;.,; for the last key value 
in all rectangles (m Xn’). Thus using the abridged 
notation from the preceding section and the 
general relation 


(10.5b) 
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FioF oF 14F isis Fi2F is Fy548 Fy4F i548 Fis.5 


we have 





(Ci) 445 = (Ci) a Fas, (a) 
(C1) 548 = (Ci) sF 54s, (b) 
>(12.1) 
(C1) 648 = (Ci) 6F eis =(Ci)s FoF oss, (c) 
(Ci) 748 = (C1) 6 Fras = (Cr) aF 6 Fris, etc. (d) 


Similarly, 





(Co) 445 = Fass, (a) 
(C2) 545 = Fosys, (b) 

(12.2) 
(C2648 = (Co) 6 Fors = FoF o+s, (c) 
(C2)748 = (Co) 6F748= FeFr4s, etc. (d) 


A summary of all coefficients for the last key 
value in rectangle (mXn’), m=3, 4, 4£n< 16 
is given in Table I. Hereafter the last key value 
will be denoted by «, in all rectangles. 


13. GENERAL RELATIONS BETWEEN THE 
COEFFICIENTS IN TABLE I 


Inspection of Table I shows three important 
relations between the coefficients for the last 
key value. 

(a) For any fixed value n’ 
(Cy) ne =(Ci)sF oF sFio- +: Fr: 

C3=C2/Fe; Cs=Cs3/Fs; 

Ci=Ci-1 F3;; C.= Fy. 


C= Cy (Cy)43 


(b) The last four coefficients C;, Cy-1, Ci—2, 
C.-; become constant for all even values of n, 
and for all odd values of n, at n212. Thus in 
rectangles m(n+ 6) for n even, 
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C;, = Fis; Ci-1 = F Fass: 

(13.2) 
C; o= F,, oF .F n+s, C; 3= F, oF oF nF ns. 
For n odd, 


Cem Fass; Cy = F 1 F 43: 
Cree = Fas Pn-1F nis; (13.3) 


C; 3= F,, 5F,, 3F 1 Fass. 


Introducing the approximations of (11.1) we 
have for n even, 


Cye=Fais; Cra = FsF sis; 
(13.4) 
C; -3 = Fe Fs.s; ts. 3= F 3 Fs... 
For n odd, 
C; 7 Fr: C; 3 = FF 743; 
(13.5) 


Cr-2= F2F743; C.-3= Fo Fr15. 


(c) For nm greater than six, assuming Fs= Fo 
= Fyo, etc. 


(Ci) nrz2= (Cia Fs. (13.6) 
raBLeE II 
u n 6- 1 oan 3 4 
vu. 2  G 0149 02128 02970 .03704 
uy 3 Ci .06748 .06818 .06878 .06931 
u 4 Ci .00082 .00152 .00213 .00265 
C. 01150 02131 .02977 .03714 
u 5 Ci .00484 .00490 .00494 .00498 
C2 .06781 .06852 .06912 .06965 
3 6 C; .00006 .00011 .00015 .00019 
C2 .00083 .00153 .00214 .00267 
C3 .01150 .02131 .02977 .03714 
U3 7 Ci .00035 .00035 .00035 .00036 
C2 .00487 .00492 .00496 .00500 
C3 .06781 .06852 .06912 .06965 
ue 8&8 CC, £00000 00001 00001 00001 
Ce .00006 .00011 .00015 .00019 
C3; .00083 .00153 .00214 .00267 
Cy .01150 .02131 .02977 .03714 
U4 9 Ci .00002 .00003 .00003 .00003 
C2 .00035 .00035 .00035 .00036 
C; .00487 .00492 .00496 .00500 
Cs .06781 .06852 .06912 .06965 
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14. MAIN TABLES AND THEIR USE 

From the expressions given in Table I, and 
the simplifications for F made in Section 11, it 
is possible to calculate the coefficients Ci, C2: - - C, 
for the 1st key value in a rectangle of arbitrary 
length mn’. Such calculations were made for 
m(n+6), m=3, 4; n=2, 3, 4, 5, 6---13, and 
5=.1, .2, .3---.8, .9, 1.0. The results are given 
in Tables II and III. 

These tables are sufficient to determine all key 
values in rectangles (m Xn’), for values of n’ 
between 2 and infinity. The procedure is as 
follows. 

First we calculate the last key value u, using 
the appropriate coefficients in the tables. We 
can then treat the rectangle as if it were shorter 
and its length m was an even integer. Next we 
find m,_, and then view the rectangle as if it 
were further shortened by 2 units. We thus can 
find 2, M3, etc. and continue the process 
until all key values are found. 

It should be observed that in calculating u,_; 
the last perimeter P;,_; must be modified by the 
addition of u,. In general, in calculating u; we 
add u;,; to P;. Thus 


u;=CyP1+CoP2+---Ci(Pituiss). (14.1) 
. 3(n+85). 

5 6 RS 8 a) 1.0 
04348  .04918 05426 .05882 06294  .06667 
06977. 07018 .07054 .07087 07116 .07143 
00312 00353 .00389  .00422 00452  .00478 
04361 .04935 .05447 .05907 06322  .06699 
00501 .00504 .00506 .00509 00511 .00513 
07012 07053. 07090 = .07123 0715207179 
00022 .00025 .00028  .00030 .00032 00034 
00313 00354 00391 .00424 .00454 00481 
04361 .04935 .05448 .05907 .06322 .06699 
00036  .00036 00036 = ,00037 ~=—-.00037 ~—- .00037 
00503 00506 .00509 00511 .00514 00515 
07012-07053. +s «07090 -~—Ss 07123-07153 ~—-.07180 
00002 .00002 00002 00002 + ~=—.00002 ~—_ .00002 
00022 00025 .00028  .00030 .00033 .00035 
00313 00354 .00391 .00424 .00454 .00481 
04361 04935 .05448  .05907 .06322 06699 
00003 00003 00003  .00003 .00003 00003 
00036 .00036 .00036 .00037 .00037 .00037 
00503 00506 .00509  .00511 .00514 00515 
07012 07053. +~=—-.07090»~—s 07123 ~—- 07153 ~—-«.07180 
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TABLE III. 4(m+84). 














u n 56— ss od oS 4 
4 
uy 2 Cy .02326 .04348 .06122 .07692 
uy 3 Ci .14667 .15000 .15294 .15556 
Ue 4 C1 .00389 .00730 .01031 .01299 
C2 .02335 .04380 .06186 .07792 
u 5 Cy .02506 .02564 .02616 .02662 
C2 .15034 .15385 .15694 .15970 
U3 6 Ci .00067 .00125 .00177 .00223 
Ce .00400 .00751 .01061 .01336 
C3 .02335 .04380 .06187 .07795 
U3 7 Ci .00430 .00440 .00449 .00457 


C2 .02579 .02639 .02692 .02740 
C3 .15045 .15396 .15706 .15982 


4 8 Ci .00011 .0002 1 .00030 .00038 
C2 .00069 .00129 .00182 .00229 
C3 .00401 00752 .01062 01337 
Cs .02335 04381 .06187 07795 


it 9 Ci .00074 .00075 .00077 .00078 
C2 00443 .00453 .00462 .00470 
C; 02581 .02642 .02695 .02742 
Cy .15045 .15396 .15706 .15982 


Us 10 C; .00002 -00004 .00005 .00007 
C2 00012 .00022 .00031 .00039 
C3 .00069 .00129 .00182 .00229 
Cs .00401 00752 .01062 01337 
Cs .02335 .04381 .06187 07795 
Us 11 1 .00013 .00013 .00013 .00013 
2 .00076 .00078 .00079 .0008 1 
"3 .00443 .00453 .00462 .00470 
4 02581 .02642 .02695 02742 
i 15045 .15396 .15706 .15982 


Us 12 Cy .00000 .00001 .00001 .00001 
C2 .00002 .00004 .00005 .00007 
C3 .00012 00022 .0003 1 -00039 
Cs .00069 .00129 .00182 .00229 
Cs .00401 .00752 .01062 01337 
Cs 02335 .04381 .06187 07795 


Ue 13 .00002 .00002 .00002 .00002 


1 2 2 

, .00013 .00013 .00013 .00014 
3 .00076 .00078 .00079 .00081 
4 .00443 .00453 .00462 .00470 
5 02581 .02642 .02695 .02742 
f .15045 .15396 .15706 .15982 


Inspection of Tables II and III shows that 
the coefficients for u, rapidly approach con- 
vergence. In rectangles (3Xn’), Table II, con- 
vergence is reached at »=9. For values of n 
greater than 9, C,, Cy-1, Cr—2, Cy_3 are the same 
as Cy, C3, C2, and Ci, respectively, at n=9 or 
n=8 depending on whether n is odd or even. 

In rectangles (4Xn’) convergence is reached 
at n=13. The coefficients for values of 6 other 
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5 6 7 8 & 1.0 





.09091 .10345 11475 .12500 13433 -14286 
.15789 . 16000 .16190 .16364 16522 .16667 


.01538 01754 .01950 .02128 .02290 .02439 
.09231 .10526 .11699 12766 .13740 .14634 


.02703 .02740 02773 .02804 02832 02857 
.16216 .16438 .16639 16822 .16990 17143 


.00264 .00301 .00334 .00365 .00393 .00418 
.01583 .01805 .02007 .02190 02357 .02510 
.09235 10532 11706 12774 .13749 .14644 


.00464 .00470 .00476 00481 00486 .00490 
.02782 .02820 .02855 .02886 02915 .02941 
.16229 -16451 .16653 .16836 .17003 17157 


.00045 .00052 .00057 .00063 .00067 .00072 
.00272 .00310 .00344 .00376 .00404 00431 
01584 .01807 .02008 .02192 02359 02513 
09235 .10532 .11706 12774 .13750 .14645 


.00080 .0008 1 .00082 .00083 .00083 .00084 
.00477 00484 .00490 .00495 .00500 00505 
.02784 .02823 02857 .02889 .02917 .02944 
.16229 16452 16653 .16836 .17004 17157 


.00008 .00009 .00010 .00011 .00012 .00012 
.00047 .00053 .00059 .00064 .00069 .00074 
.00272 .00310 .00345 .00376 .00405 .00431 
01584 .01807 .02008 .02192 .02359 02513 
.09235 10532 .11706 .12774 .13750 .14645 


00014 .00014 .00014 .00014 .00014 .00014 
.00082 .00083 .00084 .00085 .00086 .00087 
.00478 .00484 .00490 .00496 .00501 .00505 
.02784 02823 02857 .02889 02917 .02944 
.16229 -16452 .16653 .16836 .17004 17157 


00001 .00002 .00002 .00002 .00002 .00002 
.00008 .00009 .00010 -00011 .00012 .00013 
.00047 .00053 .00059 .00065 .00069 .00074 
.00272 .00310 .00345 .00376 .00405 00431 
01584 .01807 .02008 .02192 02359 .02513 
.09235 .10532 11706 12774 .13750 14645 


.00002 .00002 .00002 .00002 .00002 .00002 
.00014 .00014 .00014 .00015 .00015 .00015 
.00082 .00083 .00084 .00085 .00086 .00087 
.00478 .00484 .00490 .00496 .00501 .00505 
.02784 .02823 .02857 .02889 02917 .02944 
.16229 .16452 .16653 .16836 .17004 17157 








than those given in Tables II and III can be 
found by interpolation. 


15. AUXILIARY TABLES 


It is sometimes desirable to be able to find 
directly a value of u; +, without having to find 
first u%, Up-1°* -Ui41. This is especially useful in 
very long rectangles. This can be done by means 
of Tables IV and V, which give the coefficients 
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for any u;+u, independent of all other key 
values. Tables IV and V are for regular rec- 
tangles, with integral lengths. If the length is not 
an integer it is necessary first to find u,. 

Another advantage of Tables IV and V lies 
in the fact that all key values are given in terms 
of the one set of P’s calculated for the rectangle 
(m Xn), and no modification of these are needed. 

Inspection of Tables IV and V shows that 
when mn equals 14 the coefficients have become 
stabilized and that these coefficients may be 
used for all subsequent values of n. Specifically, 
for any integral value of m greater than 14. 

(a) The coefficients of m4, ue, and uv; are the 
same as of 4, Me, and uz for n= 14. 

(b) For » even, the coefficients C;—;, Cy—», ete. 
of um; and ug_2 are the same as of us and w,, 
respectively, at m=14. For nm odd, the corre- 


TABLE IV. (3X »). 


sponding coefficients are the same as of us; and 
us, respectively, at n=13. 

(c) The coefficients of u;, (3<i<k—2), are 
constant, and essentially the same as for us; in 
n=14. Thus, to four decimal places, odd or 
even, 


(Ui) 3m = -0722P;+ .0052(P 1+ P41) 
+ .0004(Pi2+Pis2), (a) 
(;) sen = .1768P;+.0303(P;_1+P is) 
+ .0052(Pj22+P;.2) 
+.0009(P;-3s+Piss), (b) 


since the last perimeter in (15.1) is P;+; it 
follows that u; is affected by a relatively small 
length of immediately adjacent boundaries. If 
we limit the accuracy of the coefficient to four 


TABLE V. (4X). 


or 


Ce 





u C1 C3 C4 Cc CC. 
5 wu, 07177 .00478 5 wu, .17073 .02439 
6 nu, 07179 00513 6 wm .17143 .02857 
7 wu, .07180 00515 .00034 7 My 7155 .02929 00418 
us 00515 .07214 .00481 us .02029 .17573 .02510 
8 wu; .07180 00515 .00037 8 wu, .17157 .02941 .00490 
ue 00515 .07216 .00515 us 029 .17647 = .02941 
9 wu, 07180 00515 .00037 .00002 9 ny 17157 .02943 .00503 .00072 
ue 00515 .07217 .00518 .00035 ue .02943 .17660 .03015 .00431 
us .00037 .00518 .07214 .00481 us .00503 .03015 .17588 .02513 
10 wu, .07180 .00515 .00037 .00003 10 mw, .17157 .02944 .00505 .00084 
ue 00515 .07217 .00518 .00037 ue .02944 .17662 .03028 .00505 
us 00037 .00518 .07217 .00515 uy .00505 .03028 .17662 .02944 
11 mu; .07180 .00515 .00037 .00003 .00000 11) mw, .17157 .02944 .00505 .00086 .00012 
ue 00515 .07217 .00518 .00037 .00002 ue .02944 .17662 .03030 .00517 .00074 
us .00037 .00518 .07217 .00518 .00035 us .00505 .03030 .17675 .03018 .00431 
uy .00003 .00037 .00518 .07214 .00481 us .00086- .00517 .03018 .17588 .02513 
12 wu, .07180 .00515 .00037 .00003 .00000 12 mw, .17157 .02944 .00505 .00087 .00014 
ue .00515 .07217 .00518 .00037 .00003 uo .02944 .17662 .03030 .00519 .00087 
us .00037 .00518 .07217 .00518 .00037 us .00505 .03030 .17677 .03030 .00505 
uy .00003 .00037 .00518 .07217 .00515 us 00087 .00519 .03030 .17662 .02944 
13 uy, .07180 00515 .00037 .00003 .00000 00000 13 wu, .17157 .02944 .00505 .00087 .00015 .00002 
ue 00515 .07217 .00518 .00037 .00003 .00000 us .02944 .17662 .03030 .00519 .00089 .00013 
us .00037 .00518 .07217 .00518 .00037 .00002 us .00505 .03030 .17677 .03032 .00518 .00074 
uy .00003 .00037 .00518 .07217 .00518 .00035 us .00087 .00520 .03032 .17675 .03018 .00431 
us .00000 .00003 .00037 .00518 .07214 .00481 u; .00015 .00089 .00518 .03018 .17589 .02513 
14 wu, 07180 .00515 .00037 .00003 .00000 .00000 14 wx, .17157 .02944 .00505 .00087 .00015 .00002 
ue 00515 .07217 .00518 .00037 .00003 .00000 ue .02944 (17662 .03030 .00520 .00089 .00015 
vu; .00037 .00518 .07217 .00518 .00037 . 3 ux; .00505 .03030 .17677 .03033 .00520 .00087 
us .00003 .00037 .00518 .07217 .00518 .00037 u; .00087 .00520 .03033 .17677 .03030 .00505 
us .00000 00003 .00037 .00518 .07217 .00515 u; .00015 .00089 .00520 .03030 .17662 .02944 
740 
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decimal places it is seen that in rectangles 
(4X) the length of the influencing boundary 
does not extend beyond six lattice units measured 
from point (7) in either direction. In rectangles 
(3Xm) the corresponding length is only four 
lattice units. These facts are useful in approxi- 
mate solutions. 


16. NUMERICAL EXAMPLE 


Given a rectangle 48.5 with boundary values 
as shown in Fig. 9. Required to find all key 
values. Solution : 

(a) We begin by calculating P, Pe, P3, and ps. 
The necessary Eqs. are (7.8) and (7.6). Thus 


P, =.5(4.28+3.69+2.08+3.69) 
+3.95+3.00+3.00 = 16.82, 

P.=.5(4.60+3.95+3.00+4.22) 
+4.28+3.69 = 15.86, 

P;=.5(4.914+4.28+3.69+4.64) 
+4.60+4.22 =17.58, 

ps=.5(4.60+4.22+4.91+4.64) 
+4(4.87+4.81) =47.91- 


We next find u;,=u,4. From Table III we have 
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Fic. 10. 


VOLUME 17, SEPTEMBER, 1946 


us=.00045P,+.00272P.+.01584P; 

+ .09235p4= 4.754. 
We now view the rectangle as a (4X8) and find 
U3. Thus 
uz =.00490P,+.02941P, 


Considering the remaining figure as a (4X6) we 
have 


U2 = .02857P,+.17143(P2+4u3) =3.950. 
Lastly 
uy =.16667(P)+ U2) =3.462. 
(b) Alternate procedure for u; and ue. After 


finding wu, and wu3, we take coefficients from 
Table IV and write 

















— 





uo = .02941P,+.17647P,2 

+ .02941(P3+4u3) = 3.939, 
u,=.17155P,+.02941P, 

+ .00490(P3+ U3) = 3.460. 


The exact values of 44, 43, M2, and “; are 4.7536, 
4.3820, 3.9512, 3.4657, which differ but slightly 
from those given by our method. From the key 
values all remaining values can be calculated by 
the Liebmann formula. 


17. COMPOSITE RECTANGULAR AREAS 


(a) Angles. The proposed method and tables 
are applicable to composite rectangular areas. 
Consider the angle of Fig. 10. This configuration 
may be viewed as the composite of two rec- 
tangles (mXn) and (m Xm). The key values 
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Uy, Ug ++, in rectangle (m Xn) could be found 
by means of the coefficients in the tables if u»’ 
were known. Similarly 1“, u2’---u,’ in the rec- 
tangle (mXn,) could be calculated if uz. were 
known. The necessary key values uz and u,’ can 
be found from the following two equations: 


ue = Cyi(P\ +2’) +C2(P2+.5u2’) 


+C3P3+---C,Px, (Db 097.4) 
us’ = Cy’ (Py +2) +Ce' (Po! +.5u2) , 
+C;'P;'+ ee Cy’ Py’, (b) | 


in which the coefficients C and C’ are taken 

directly from the tables, and most perimeters are 

calculated from the basic expression 
P=.5 > known a’s+)> known c’s, (17.2) 


using only the boundary values on the angle. 
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We denote the results from (17.1) by 


u2=f(C, C’, P, P’), (a) (17.3) 
us’ =f,(C, C’, P, P’). (b) on 


(b) T-figures. A configuration such as shown 
in Fig. 11 may be viewed as the composite of 
_ two overlapping angles. Observing that P; in 
angle A equals P, in the T-figure plus u.’, and 
that P,’ in angle B equals P,’ in the T-figure 
plus “2 we have by (17.3) 


ur=f(C, P, C", P”’, us’), 


Ue’ =f,(C’, | . . Us). (17.4) 


Upon solution the problem is reduced to angles. 

An alternate procedure is to view Fig. 11 as 
consisting of three overlapping rectangles of 
length m, m,, and me. The resulting equations are : 
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Uo=Ci(Pitue’+u2"’) ) 
+C2(P2+.5u2!)+C3P3+---C.P:, (a) 
us’ = Cy'(Py'+u2+u2’’) 
+ C2! (Po! +.5u2’’) 
+C3'P3'+---C,’P,’, (b) 
uo” = Ci" (Py +u2+U2’) 
+ C2!" (Po! + .5u2+.5u2’) 
4+-C," Pf’ + ee CP Pia. (c) 


Upon solution the figure is reduced to rectangles. 

(c) Channels. A channel, Fig. 12, may be 
viewed as the composite of two overlapping 
angles A and B. We observe that 


(Pw) 4a™ (Par)channel +Y; (a) | 17 6) 
( Py) in (Py) channel +X. (b) s 


Hence by (17.3) 


x= f(C,C’, P,P’, y), (a) 
y=fi(C,C’, P, P’,x). (b) 


Upon solution we obtain x and y and the problem 
is reduced to angles. When the length »>3m 
the value y may be neglected in the calculation 
of x, and x may be neglected in the calculation 
of y. Equations (17.7) then reduce to (17.3). A 
figure having a Z-shape can be treated in the 
same manner. 

(d) H or I shapes may be viewed as over- 
lapping channels. The method is also applicable 
to rectangles of variable widths, Fig. 13. It 
should, however, be added that in complicated 
areas the method of block iteration may be 
shorter.’ 


>(17.5) 





}«a7.7) 
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Erratum 
Correction of Diffraction Amplitudes for Lorentz and Polarization Factors 


[J. App. Phys. 17, 285 (April 1946)] 


M. J. BUERGER AND GILBERT E, KLEIN 


HE values of the functions in Tables III and VIII are each misplaced in 

the region 90-179.9° so that they are found at a value of the given 
argument 0.1° lower than the correct argument. The tables may be utilized 
in this region by subtracting 0.1° from the desired value of T and looking up 
the function for this altered (false) argument. For example, if one desires the 


function for 147.3°, it should be looked up under the altered (false) argument 
147.2°. 





The Thermoluminescence and Conductivity of Phosphors 


RoBerT C. HERMAN AND CHARLES F. MEYER 
Applied Physics Laboratory, Johns Hopkins University, Silver Spring, Maryland 
(Received May 3, 1946) 


A simple theoretical treatment is given for the emission of light and associated variations 
of dark current that occur on warming zinc sulphide and zinc silicate type phosphors pre- 
viously irradiated with ultraviolet light at low temperatures. Reasonably good agreement is 
found between the calculated and observed glow curve for Willemite. The same type of treat- 





ment is applied to the case of infra-red illumination of phosphors at low temperatures. 


1. INTRODUCTION 


HE phenomena of thermoluminescence and 

the associated increase in conductivity in 
certain phosphors such as zinc sulphide and zinc 
silicate are sufficiently interesting to warrant 
further study. Recently there has arisen a great 
deal of interest in the storage and release of 
energy in solids, particularly with respect to infra- 
red-sensitive phosphors. Simple models based on 
the energy-band theory of solids have been 
shown to be in agreement with experimental data 
by various investigators.'~* It is now generally 
accepted that a phosphor is an insulating crystal 
which can be made luminescent by the addition 


1F, Seitz and R. P. Johnson, J. App. Phys. 8, 84, 186, 
246 (1937). 


2N. F. Mott, Trans. Faraday Soc. 34, 822 (1938). 
°R. P. Johnson, J. Opt. Soc. Am. 29, 387 (1939). 
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of a small concentration of impurity or activator 
atoms. Figure 1 shows the spectrum of energy 
states of electrons in the crystal. Between the 
uppermost filled band, A, and the next highest 
band, B, the conduction band which is normally 
empty, there are localized electron states, such 
as C, D, and E. The narrow energy states C 
and D which lie in the forbidden region are the 
impurity states associated with the activator 
atoms whereas the states E represent the so- 
called trapping states which are associated with 
lattice imperfections such as strains, cracks, 
surfaces, etc. The trapping states, E, are popu- 
lated with electrons to an extent which depends 
on the previous history of the phosphor. Photo- 
conductivity arises when an electron is raised 
into the conduction band from states A, C, or E 
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Fic. 1. Energy states of electron in phosphor. 





by the absorption of a quantum of light. How- 
ever, no photo-conductivity is observed when an 
electron is raised to an excited state at D. 

In an excited phosphor, electrons return to the 
ground state by dropping from the conduction 
band to a state at D and then may make a 
radiative transition from D to C. Electrons in 
the conduction band may also be captured in 
the trapping states E. They can be released 
either thermally or by the absorption of radiation 
and thus return to the conduction band. The 
phenomena of phosphorescence and thermo- 
luminescence are caused by the storage of elec- 
trons in excited states of the activator atoms and 
in the trapping states. Recently, Randall and 
Wilkins‘ have discussed the phenomenon of 
thermoluminescence or glow and presented a 
series of significant experiments. The storage of 
light at low temperatures which can be recovered 
by warming is well known and has been dis- 
cussed frequently. The glow curve for zinc 
silicate obtained by warming a phosphor or 
crystal previously illuminated at liquid air tem- 
peratures with ultraviolet light is accompanied 
by a large dark current peak at about the same 
temperature as the glow peak.® This observation 
demonstrates that electrons thermally released 
from trapping states pass through the conduc- 
tion band during the glow process. This con- 


‘J. T. Randall and M. H. F. Wilkins, Proc. Roy. Soc. 
A184, 366 (1945). 


®°R. C, Herman and R. Hofstadter, Phys. Rev. 57, 936 
(1940). 
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clusion is supported by the fact that a glow is 
produced in phosphorescent zinc silicate by the 
application of an electric field.6 As previously 
suggested, it would prove interesting to observe 
the variations of dark current in zinc-silicate 
type phosphors during the glow experiment.® 
Also, dark current measurements would permit 
a study of trapping state distributions in non- 
fluorescent solids.‘ 

It is the purpose of the present paper to discuss 
theoretically both the emission of light and 
variations in conductivity that occur in phos- 
phors of the zine silicate and zinc sulphide 
varieties. This has been done with a simplified 
model in which the only energy states considered 
are those in the uppermost filled band, the next 
highest normally empty conduction band and in 
the trapping state distribution. For the sake of 
simplicity, the analysis which follows deals with 
discrete sets of trapping states. This may be 
extended to will be 


various distributions as 


mentioned later. 


2. LUMINESCENCE AND CONDUCTIVITY 
DURING WARMING 


Let N; be the number of trapping states per 
cm’ at an energy ¢; below the conduction band, 
no the number of positive holes per cm* in the 
uppermost filled band, » the number of electrons 
per cm* in the conduction band, and n; the 
number of electrons per cm’ in the trapping 
states in de; at e;. If a phosphor of the type under 


discussion is irradiated with light of the proper 


wave-length at a sufficiently low temperature, 
the trapping states will become populated and 
remain so for considerable lengths of time so 
long as exp’ (—e;/RT) is very small. With regard 
to the warming experiment, the initial condition 
of the crystal or phosphor chosen is that the 
number of electrons trapped in the states in de; 
at an energy e; below the conduction band is (m;)o, 
n=0 and (m)o the number of holes, originally 
equal to the number of trapped electrons. In 
general, m»=n-+n;. If the crystal is allowed to 
warm up the light emission and variations in 
dark current may be found from the following 
set of equations: 


6 J. Frenkel, Tech. Phys. U.S.S.R. 5, 685 (1938). 
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dny/dt= —Angon, (la) 


dn 
= —Anygnt+)> s; exp (—e;/kT)n; 
at 
—> Ki{(ni)o—n; jn, (1b) 
dn; 
= —s, exp (—e,/kT)n; 
dt 


+K {(nio—n; jn, (1c) 


where A is the recombination coefficient and is 
equal to ol; o the cross section for capture by 
an impurity center of an electron in the conduc- 
tion band and subsequent radiation; V the ve- 
locity of the conduction electron; s; the frequency 
factor in sec.—' or the number of times per second 
an electron strikes the potential barrier in a trap; 
and K; is similar to A but refers to the various 
sets of trapping states. De Groot’ has solved 
graphically a similar set of equations in his 
studies of the fluorescence and phosphorescence 
of zinc sulphide phosphors. Following Randall 
and Wilkins* who-point out that the effect of 
retrapping is not very large, we have solved the 
equations for the case K;=0 or: 


dny/dt= —Anon, (2a) 
dn . 
~=—Anomt+)> s;exp (—e:/kT)n;, (2b) 
dt 
dn; 
—-=—s, exp (—e;/kT)n,. (2c) 
dt 


Equation (2c), which expresses the rate at which 
trapped electrons leave the states in de; at €;, can 
be integrated yielding: 


t 
nN, =(nj)o exp | -s.f exp (—e/ATa (3) 
0 


where (;)9 is the number of electrons originally 
trapped in de; at ¢;. If a constant rate of warming 
is used so that 7’ =at+T7>5, then dT =adt and 
n,=(N;j)o 

T 


Xexp ; —(s, a) f exp (—e;/RT)dT}. (4) 
To 





The general case involving electrons trapped 


in complicated trapping distributions can be 
7W. de Groot, Physica 6, 275 (1939). 
8]. T. Randall and M. H. F. Wilkins, Proc. Roy. Soc. 
A184, 390 (1945). 
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handled by a numerical integration of Eq. (4) 
over all trapping states. 

In the case where there are sets of states at 
discrete energies the problem is much simpler 
and Eq. (4) becomes for a single trap energy: 


n,=N, exp ie ak) g*s ds, (3) 
z0 


where 


x=e/kT, Xy9=e/kTo, 


and it is assumed that N;, traps at the energy « 
are filled with electrons when T=T7>. Equation 
(5) can be rewritten in the form 


n\= N; exp { wie (se /ak) 


X[x-'e*+ Ei(—x)—fo]}, (6) 


— k1(-—x) -{ u—'e-“du 


a 


where 


and 
fo =x! exp (—Xo9) + Fi( —Xo). 


Substitution of 2; from Eq. (6) into (2a) yields 
an equation that can be expressed in terms of 
known functions, but which requires numerical 
solution. It was found convenient to solve the 
differential equation numerically in the form: 


dn AN 7 
=" G-9), (7) 


dx ak x? 
where 


n=no/Ni, €=m,/N, and §t=n/N,=n-¢. 


The solution of these equations yields both the 
luminescent intensity, 7, and the conduction 
current, 7, in an applied electric field since J « né 
and i &, 

Randall and Wilkins* have shown that for 
many of the phosphors the glow curves can be 
interpreted on the basis of a reasonably broad 
distribution of trapping states. For example, the 
broad glow curves for zinc sulphide phosphors 
originate in this manner. However, the three 
glow-curve peaks for Willemite+1 percent Mn 
seemed sufficiently sharp to be accounted for by 
a set of three discrete trap energies. 

The light intensity and current have been 
computed for this case by solving Eq. (7) 
numerically using the value of the constants 
shown in Table I. The solutions to these equa- 
tions are shown in Fig. 2 where n, £ and ¢ are 
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Fic. 2. Curves showing temperature dependence of 
luminescence and current during linear warming of a 
Willemite phosphor in an electric field, for the case of 
peak C (see Fig. 3). 


plotted against x for the largest peak labeled C 
in Fig. 3. It is to be noted that the experimental 
fit seems to be good, particularly in that the 
long tail is predicted. For long times when the 
number of trapped electrons is small (or .¢—0) 
Eq. (7) reduces to a form having a solution of 
the type 

n=&=([x/(1+cx) ]. (8) 


Thus, for long times bimolecular decay applies. 
In the event that positive holes are originally 
present the decay would be more rapid and in 
fact exponential if there are a large number of 
holes present originally. 

The conduction current is proportional to é 
and has not been plotted on a temperature scale 
since complete experimental data is not available 
for comparison. However, Herman and Hof- 
stadter® found a sudden increase in dark current 
at about —60°C which agrees with the predicted 
‘curve shown in Fig. 2. 

Since the areas under the light emission curves 
are proportional to the number of trapped 
electrons, the relative number of trapping states 
at the three energies could be obtained. Since 
the peaks A and B overlap, it would be necessary 
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TABLE I. 

Constants Peak A Peak B Peak C 
A(cm*/sec.) 1078 1078 10-3 
s(sec.~!) 10° 10° 10° 
e(ev) 0.223 0.291 0.500 
Ni(cm~) 1.4 10" ~5X 10" 1x 10'8 
a(deg./sec.)* 5 5 3 





* The experimental data of Randall and Wilkins (reference 4) shows 
a non-uniform rate of warming for this case. 








to solve Eq. (7) in the form 


dn 


wl 1) 
ea eae — 1), 
aT _ an ( 


where N¢(T) is the number of electrons trapped 
in both sets of levels. However, for the sake of 
simplicity each set of trapping states was 
handled separately and the result is shown in 
Fig. 3. Since overlapping light emission and cur- 
rent peaks are more difficult to analyze, it would 
be desirable to have a slow linear increase in tem- 
perature since this will lead to better resolution. 
This has been done by Williams and Eyring.® 


3. LUMINESCENCE AND CONDUCTIVITY DURING 
INFRA-RED ILLUMINATION 


Some time ago Gisolf'® performed a series of 
qualitative experiments with ZnS_ phosphors 
which showed that electrons in trapping states 
absorb infra-red radiation. The electrons can be 
freed by infra-red illumination and since changes 
in conductivity occur it is evident that the 
electrons pass through the conduction band. On 
first illuminating the phosphor the photo-current 
increases and then decreases as a consequence of 
recombination of electrons with positive holes. 
A more quantitative study by Reimann" shows 
that the resistance of a ZnS phosphor containing 
trapped electrons increases linearly with time 
when the phosphor is irradiated with infra-red 
light. The simple theory for sufficiently long time 
after beginning the illumination has been given 
by Reimann and others." '% 





*F. E. Williams and H. Eyring, ‘Mechanism of lumines- 
cence of solids,” Bull. Am. Phys. Soc. 21, 17 (A) (Jan. 24, 
1946). We wish to thank Professor H. Eyring and Mr. 
F. E. Williams for interesting discussions. 

1 J. H. Gisolf, Ned. T. Natuurk 5, 289 (1938). 

11 A. L. Reimann, Nature 140, 501 (1937). 

2N.F. Mott and R. W. Gurney, Electronic Processes in 
Ionic Crystals (Oxford University Press, New York, 1940). 

13F, Seitz, The Modern Theory of Solids (McGraw-Hill 
Book Company, Inc., New York, 1940). 
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Fic. 3. Comparison of theoretical glow curve with 
experimental results of Randall and Wilkins (reference 4) 
for Willemite +1 percent Mn. 


It is evident that the resistance of such a 
phosphor must decrease and then increase in a 
manner which depends upon the infra-red light 
intensity, the cross section of a trapped electron 
for an infra-red photon, the total number of 
trapped electrons, and the probability that a 
conduction electron will recombine with a posi- 
tive hole. Since a general theory involving com- 
plicated trap distributions leads to tedious solu- 
tions, attention has been confined to the case of 
a single set of trapping states all at the same 
energy below the conduction band. 

The equations are similar to those employed 
in the previous case except that the temperature 
is kept constant and the rate of release of elec- 
trons from trapping states is given by 


dn, Io 
——~ 2S arene. (9) 


dt hv 


where Jo is the intensity of the infra-red light, 
o the cross section of a trapped electron for an 
infra-red photon, and m, the number of trapped 
electrons. For constant intensity of illumination, 
To, Eq. (9) yields 


n)i= N, exp (— Col), 


Co= Too. hv, 


(10) 


where 


and N, is the number of electrons originally 
trapped. Substitution of Eq. (10) into (2a) 
yields 


dno 

= —Anlno— MN, exp (— Cot) ]. (11) 
The solution of Eq. (11) is 

1 AN, 

-=1———[Ei(—y)—Ei(—yo)]e", (12) 





n Co 
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Fic. 4. Time dependence of luminescence and current 
for a phosphor in an electric field, that is illuminated with 
infra-red light at constant temperature. 









































where 
y=(AN,/Co) exp (—Cot), you ANi/Co 


and 


n=No/ N,. 


When ¢ is sufficiently large §=(n/N;) is equal to 
n and the resistance which is proportional to 
(1/£) increases linearly with the time. 

The solutions for the glow intensity and the 
current are shown in Fig. 4 for the case of 
AN,=1 sec.! and Cy=0.1 sec.~!. The equations 
have also been solved for various values of AN, 
and C, and the resistance computed as a function 
of time. These results are shown in Fig. 5. It is 
apparent that if the infra-red intensity or the 
cross section is smaller, the resistance minimum 
is broader. On the other hand, if the total num- 
ber of electrons originally trapped or the recom- 
bination coefficient is smaller, then the resistance 
begins to rise linearly at an earlier time. Since the 
resistance in the dark is finite, it is necessary 
to add a constant term to the conductivity 
calculated. However, if the peak current is large 
compared to the dark current, deviations from 
linearity will occur only at very long times. 
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Fic. 5. Variation of resistance «1/£ with time, for a 
phosphor in an electric field during infra-red illumination 
and the dependence on phosphor characteristics as well as 
light intensity. 


4. APPENDIX 


The exponential integral 


m ” e u“ 
- bi(—x)= f —du 
u 


£ 


employed in the calculation of the glow intensity 
curves, etc., has been tabulated in various 
_ references'*"'® for OSx=15 and 20=x=50. We 
i Tables of Sine, Cosine, and Exponential Integrals 
(Federal Works Agency, WPA, 1940), Vol. 1 and 2. 

'® British Assoc, Adv. Sci. Math. Tables (Cambridge 
University Press, New York, 1931), Vol. 1. 

‘6 T, Akahira, Sci. Pap. Inst. Phys. Chem. Res. Tokio, 
181 (June, 1929), Table No. 3. 
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TABLE II. Exponential integral — Ei(—x)= f.°(e~"/u)du. 


x —Ei( —x) As 
15.0 1.91863 x 107% 0.02178 1075 
15.1 1.72520 1078 0.01956 x 10° 
15.2 1.55133 107% 0.01756 10°* 
15.3 1.39502 K 10° 0.01581 «107 
15.4 1.25452 X10°5 0.01419 10 3 
15.5 1.12821 1075 0.01276 10 8 
15.6 1.01466 x 10°* 0.01146 10° 
15.7 9.12573 X10" 0.01030 x 10-5 
15.8 8.20784 x 10°° 0.09013 x 10°* 
15.9 7.38008 x 10~* 0.08818 x 107° 
16.0 6.64050 x 10~° 0.07229 x 10-° 
16.1 5.97321 10~° 0.06727 x 10~° 
16.2 5.37319 x 1079 0.06042 x 10~° 
16.3 4.83359 x 107% 0.05434 x 10° 
16.4 4.34833 x 107° 0.04886 x 107° 
16.5 3.91193 x 107° 0.04390 x 1079 
16.6 3.51943 x 107° 0.03949 x 10-9 
16.7 3.16642 x 10~° 0.03550 x 107° 
16.8 2.84891 x 107° 0.03191 x 107° 
16.9 2.56331 10-9 0.02872 x 10-9 
17.0 2.30643 X 107° 0.02580 x 107° 
17.1 2.07535 X 107° 0.02321 x 10-° 
17.2 1.86748 x 10~° 0.02088 x 10-° 
17.3 1.68049 x 10° 0.01876 x 10° 
17.4 1.51226 10-° 0.01688 x 10~° 
17.5 1.36091 x 107° 0.01519 x 107° 
17.6 1.22475 « 107° 0.01366 x 10-° 
17.7 1.10225 10~° 0.01227 xX 107° 
17.8 9.92022 x 10-" 0.01106 x 10~° 
17.9 8.92850 x 10-” 0.10031 107" 
18.0 8.03609 x 10~" 0.08843 x 107° 
18.1 7.23311 10-" 0.08041 x 10~"” 
18.2 6.51054 x 10°" 0.07233 x 10~” 
18.3 5.86030 x 107" 0.06509 x 10-" 
18.4 5.27515 x 10-” 0.05859 x 10-"° 
18.5 4.74859 x 107° 0.05263 x 10-” 
18.6 4.27466 x 10-” 0.04741 x 10-" 
18.7 3.84814 x 107” 0.04267 x 10-"” 
18.8 3.46429 x 10-” 0.03835 X 10~-” 
18. 3.11879 x 10-” 0.03454 x 107° 
19.0 2.80783 x 10-” 0.03042 x 10-” 
19.1 2.52729 X 10-” 0.02924 x 10-” 
19.2 2.27599 x 10-” 0.02453 x 107° 
19.3 2.04922 x 10-” 0.02263 x 10-” 
19.4 1.84508 x 10-” 0.02037 x 10-” 
19.5 1.66131 x 10-” 0.01835 x 10-” 
19.6 1.49589 x 107° 0.01651 x 10-" 
19.7 1.34698 x 107° 0.01484 x 10-” 
19.8 1.21291 X10" 0.01337 x 10-” 
19.9 1.09221 x10-” 0.01204 x 10~" 
20.0 9.83555 K 10-4 0.01083 x 10-” 





have been unable to find tabulated values of 
—FKi(—x) in the range 15<x<20. These have 
been calculated using the asymptotic expansion 
ae oe | ae ae 
~ HM ~2) 06°94 --— 4 — ——  — « « 
x x 2 x 3 x 4 x 
and are given in Table II together with second 
differences, As, in the hope that they may be 
found useful by others. 
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Oscillation Conditions in Single Tuned Amplifiers 


WILLIAM R. Faust AND HucGo M. BEck 
Naval Research Laboratory, Washington, D. C. 


(Received May 31, 1946) 


An application of the calculus of finite differences has been made in obtaining the complete 
expression for the voltage gain of an n-stage amplifier having identical grid-to-plate impedances 
and plate loads when driven by a generator of any given internal impedance. An application 
of these results was made in determining the conditions of oscillation in a multi-stage amplifier 
using single-tuned circuits. It was found that there is a certain minimum grid-to-plate capaci- 
tance required to cause oscillation. If the actual grid-to-plate capacitances are below this 
value, the system will not oscillate. There also exists a region of stable gain, zero to 2", wherein 
it is impossible to make the amplifier oscillate even with grid-to-plate capacitance arbitrarily 


large. 


1. INTRODUCTION 


HE calculus of finite differences has been 

applied to cascaded systems of vacuum 
tube amplifiers,| but apparently the full power 
and usefulness of the method has not been 
realized. It is obvious that a system of cascaded 
amplifiers lends itself quite readily to the methods 
of finite differences. The equivalent circuit of 
multitube amplifiers resembles, to a very great 
extent, an artificial transmission line made up of 
Pi or T sections with the exception that the 
amplifier contains internal current generators. 
Nevertheless, the difference equations for the 
two systems are quite similar in form and become 
identical when the g» (mutual conductance) of 
the tubes vanishes. 

The purpose of this paper is twofold: First, 
that of obtaining a complete expression for the 
voltage gain of a multi-stage amplifier. Such an 
expression describes the effect of the feedback 
between the grid and plates of the various stages. 
Second, that of applying the general expression 
for voltage gain to specific cases and using it to 
compute the conditions necessary for sustained 
oscillations to exist in a multi-stage, single-tuned 
amplifier. 


2. DEVELOPMENT OF FUNDAMENTAL RELATIONS 


The circuit of Fig. 1 is a cascaded multi-stage 
amplifier driven by a generator of e.m.f. e;, and 
internal impedance Z, with identical plate load 
impedances Z,’. The circuit of Fig. 2 is an 

' A. C. Bartlett, The Theory of Electrical Artificial Line 


and Filters (John Wiley and Sons, Inc., New York, 1931), 
pp. 143-150. 
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equivalent circuit of Fig. 1, wherein Z,; represents 
the grid to plate impedance which is assumed to 
be perfectly general at the present time. For 
simplicity the circuit of Fig. 2 will be shown as 
in Fig. 2a, where Z» is the parallel combination 
of R, and Zy’, i.e., 


Z2=R,Z2' /(Z2' + Rp). 
The voltage relationship between the k and 
k+1 stages is 
C41 —-Zik —e. = 0, (1) 
and the currents at the plate of the kth stage are 
related by 
ty, — 1-1 — €/Z2— Zmex—-1 = 90. (2) 
Equations (1) and (2) are the difference equations 
of the amplifier for which a solution satisfying 
the boundary conditions must be found. To find 
a solution, it is assumed that 
e,. = U(exp yik) + Viexp —vyek), (3) 
1, = W(exp yik) +X (exp —72k), (4) 
where UL’, ’, W, and X are complex constants and 
y, and y2 are the complex propagation constants 
of the system. There are not four independent 

















Fic. 1. Radiofrequency circuit of m-stage amplifier. 
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Fic. 2. Equivalent circuits. (a) Equivalent circuit of Fig. 1. 
(b) Simplified equivalent where Z2= R,Z2'/(Z2'+R>). 


constants U, V, W, X but only two, as a relation- 
ship between U and W, and V and X can be 
found. By substituting the assumed solutions (3) 
and (4) into Eq. (1), it follows that 


(exp yk) U(—1+ exp y1) —Z,W] 
+ (exp — yk) V(—1+exp —y2)-—Z:X ]=0. (5) 


Since this expression must hold for all values of k, 
it is necessary that the coefficients of exp y:k and 
exp y2k vanish independently, i.e., 


W= U[(exp y1)—1]/Z,, (6) 
X = V[(exp — 2) —1]/Z:. (7) 
Equation (4) then becomes 


i, = U(—1+ exp y1) 
X (exp yik)/Zi+ V(—1+exp —vy2) 
X(exp —ye2k)/Z;. (8) 


To determine the “propagation constants”’ y; 
- and 2 so that the assumed solutions will satisfy 
the difference equations, substitute the solutions 
(3) and (8) into Eq. (2), which then becomes 


U{(—1+exp y: ]Lexp y:k]/Z1 
—[—1+exp y: ][exp yi1(k-1)]/Z; 
+ —[exp vik ]/Z2—gmLexp yi(k —1)]} 
+V{l—1+exp —y2][exp —y2k]/Z; 
+—[—1+exp —y2]Lexp —y2(k—-1)]/Z1 
—[exp —y2k]/Z2 
—mLexp —y2(k—1)]}=0. (9) 
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Because this must be satisfied for every value of 
k, it is necessary that the coefficients of U and |’ 
vanish; and since UV’ and V are assumed not to 
equal zero, it follows that 


exp y1=1+2Z,1/22Z2 


+[Z1/Z2+(Z1/2Z2)?+gmZi]', (10) 
and 
exp y2=1+2Z,/2Z, 
+[Z, Z2+(Zi 2Z2)* + gmZ1 |*/(1—gmZ1). (11) 


The exact choice of the signs in (10) and (11) will 
be determined later. It should be pointed out, 
however, that (exp yi)#(exp —y2) as such a 
choice of signs would lead to a solution not 
consistent with the physical facts. It has been 
shown that Eqs. (3) and (8) are solutions of the 
difference equations where y; and ye are as 
defined above and U and V determined by the 
boundary conditions. 


3. BOUNDARY CONDITIONS 


The boundary conditions are very similar to 
those for a transmission line and should require 
but little discussion. It will be assumed that the 
amplifier has m stages, so the corresponding 
boundary condition is simply 


i, =0, (12) 
from which it follows that 
V=—U[exp n(yi+72) ] 
x[—1+exp y1]/(—1+exp —y2). (13) 


Since the amplifier is driven by a generator of 
e.m.f. ej, and internal impedance Z,, the corre- 
sponding boundary condition is 


Co =Cin tZoto. (14) 


By the use of Eqs. (1), (8), and (13) Eq. (14) 
becomes 


ein = U\1—[Lexp u(y +72) ] 
X{—1+exp y:]/(—1+exp — 2) 
+—Z,(—1+exp 71) 
x1 —exp n(yitve) |/Zi}. (15) 
Consequently, the expression for the voltage gain 
between the generator and the kth stage becomes 
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e; Lexp yk ]—Lexp n(yi+v72) ][—1+exp yi [exp —y2k]/[—1+exp —72] 


=——— 3 . (16) 
ein 1—{Lexp n(yi+v2)][—1+exp y:]/[—1+exp —y2]} 





Z 
ae rs Se 1+exp y1][1—exp n(yi+72) ] 


1 


Usually the over-all gain of the system is of interest rather than the gain up to a certain stage, so 
the expression for e,/e;, will be considered. This is obtained from (16) simply by setting k=n, 


Cn (exp —Y2)— (exp +71) . 
- = ——_ . (17) 
ein [exp —nyi][—1+exp —y2]—Lexp ny2][—1+exp 71] 





Zo 
“— 1+exp yi ]L—1+exp —y2]L(exp —n71) — (exp ny2) J 


1 


4. SPECIAL CASES 


Equation (17), as it stands, is of little practical value; and the exponentials must be replaced by 
their values in terms of the parameters of the system. For simplicity in the manipulation of the 
expression, it will be convenient to set 


A=! +Z;/2Z:2, 


and 
B=[Z, Z2+(Z1/2Z2)*+2mZ1 |}, 
so that 
expyi=A-—B or exp —y:=(A+B)/(1—Zign), 
and 


exp y2=(A—B)/(1—gmZ1). 


The relationships may be simplified by setting Z, = Z2. This is the usual case and should not detract 
from the utility of the results. Substitution of these values in (17) produces 











TABLE I, 
No. of stages, 1 Gain 

0 1 
1 Z2—8mZ 22) . 

2Z2+ZitgmZ2 
' (1 gna Ps* 
7 (Z14+-3Z2+22mZ2?)(Z1+Z2) 
3 (1 = gmZ1)°Z2° 


Z8+3Z2Z2(2+gmZ2)+Z1Z2(104+8gmZ 2+ Sm?Z2)+Z2(44+3gmZ2) 
( 1 s. LmZ1 YZ -4 


— Ua gmei PZ 
ZP+5Z1'Z2(2 + BmZ2) +22 3Z27(18+ 16gmZ2+32m2Z22) +ZPZF(gmiZ 2+ 50+ 1 Bgm2Z 2+ 038mzZ2) 
+5Z:Z24(7+8gmZ2+2gm?Z 22) +Z29(6+5emZ2) 
‘ ae wee (1—gmZ1)8Z28 
(ZB +7ZZ2+4gmZ2Z 2+ 14Z:Z 2+ 13g¢mZ1Z 2+ 2Em2Z1Z2+Z+OgmZ2') 
X (Zi +52°Z2+2gmZ 1°22? +621:22? + 3gmZ1Z2°+Z:) 


wn 
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Cn 2B(1—gmZ1)" 





2, 
(A+B)"(A+B—1)—(A—B)"(A —B-1)—(A—B~=1)(A+B-—1)[(A +B)"—(A —B)"}- 


which may be written in the form 
Cn Cin 
By expanding Eq. (19), it is easy to obtain 


€./€.n™ 
n(n—1)A”"-*B? 
(n-+1)| An —— + 


3! 


+ enZ| n+ ' 
3! 


~ (A+B)"((A+B)+£nZ2]—(A—B)"[(A—B) +gnZ2] 


(1 —2mZ1 - 


———— —, (18) 
“i 
2B(1—gmZ1)" 
(19) 
: (20) 
n(n—1)(n—2)(n—3)A"—*B! 
See...) 
5! 
n(n—1)(n—2)A"-*B* n(n—1)(n—2)(n—3)(n—4)A”>B! 
a eee. 
5! 
Z2=R/(1+JQP), (22) 


Equation (20) is the general expression for the 
gain of a cascaded n-stage amplifier and, as such, 
is much too general to deduce any results. In 
Table |, Eq. (20) has been evaluated for various 
values of n; the appropriate values of A and B as 
given previously have also been included. 

That these equations are correct in the limit 
may be shown by letting Z,; become very large; 
the general form of the expressions in Table | 
then becomes (— gmZ2)", which is the proper form 
for the gain of an n-stage amplifier. 


5. CONDITIONS FOR OSCILLATION 


It is now desired to determine what conditions 
must exist in order to produce sustained oscilla- 
tion. Obviously the system under consideration is 
oscillating when a finite output e, #0 exists for 
ein =0. In order that the above conditions obtain, 
it is necessary that the denominators of the 
various expressions in Table | vanish. To apply 
the above results to a practical case, Z, will be 
taken as the plate to grid impedance of a vacuum 
tube which is mainly capacitive, i.e., 


Z1=1/juC. (21) 


The effective plate impedance Z, will be taken 
as that of a single tuned circuit consisting of an 
inductance Lo, capacitance Co, and resistance, all 
in parallel; hence? 


2S. A. Schelkunoff, Electromagnetic Waves (D. Van 
Nostrand Company, Inc., New York, 1943), pp. 119-120. 
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where 
QV= R( Co Lo)}, 
F= (\w/wo— wo/w), 


wo —— . : 
5>— = Resonant frequency of single-tuned circuit 


== | 2x(LoCo)}, 


Ww _— ° . . 
,~ = The frequency under consideration, and 
- : 


R= Dynamic resistance 


= R’R,/(Rp+R’) ~R’ (for pentodes). 


One Stage 


From Table |, it is seen that the conditions of 
oscillation are to be determined from 


2Z2+ZitEmnZ2 =0. (23) 


By substituting the values of Z; and Zs from 
Eqs. (21) and (22) and separating the real and 
imaginary parts, there is obtained 


(QF)?+2(RwC)(QF) —1=0 (Real part), (24) 
and 
2(QF) +2(RwC) 
—gnR(RwC)=0. (Imaginary part). (25) 
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A simultaneous solution of these two expressions 
results in 


QF = +[(2+8mR)/(gmR—2) }}, (26) 


and 


C= ¥2/wR[(gmR)?—4]}. (27) 


The negative sign must be chosen in Eq. (26) so 
that C will be positive in Eq. (27). The frequency 
of oscillation w is determined from 


OF = (w/w — wo/w)Q = (w? — wo?) Q/ ww 


= —[(2+gmR)/(gmR—2)]?, (28) 


which gives 
_ Wo (gmR+2)7} 
“a oa 
7 pene 


" 
+40] ; (29) 
(gnR—2) 


The positive sign of the second term in (29) must 
be chosen so as to produce a frequency of 
oscillation that has physical meaning. After w is 
eliminated from Eqs. (27) and (29) there results, 


40 





~ Reool [(gmR +2)? +402(gn?R?—4) } — (gnR+2)} 


It is interesting to note that if gnR<2, the 
system will never oscillate, as C becomes indefi- 
nitely large when g,,.R=2, and is imaginary if 
gn <2. For a high stage gain, these results differ 
but little from previous results.* * For an amplifier 
of low gain (g,R of the order of 2) and low Q, 
there is a wide difference; and these results indi- 
cate that the system will be intrinsically stable. 
This situation is illustrated in the curves of 
Fig. 3 in which Eqs. (24) and (25) are plotted 
with (QF) as the independent variable and (g,,R) 
a parameter. Equation (24) is a hyperbola and 
(25) is a straight line whose slope is dependent 
upon the parameter g,,R. The intersection of 
these two curves represents the simultaneous 
solution of the real and imaginary parts of the 
defining equation. In order to obtain positive 
values of (RwC), the intersection of the two 
curves must occur in the second quadrant. Since 
the asymptote to the hyperbola corresponds to 
the straight line with g,,R = 2, it follows that there 
can be no real simultaneous solution of the two 
equations between g,R=0 and g,,R=2. Conse- 
quently, the system will be inherently stable 
within these limits. 


Two Stages 
The defining equation in this case is 


(Z1+Z2)(Zi1+3Z2+ 2gmZ2*) =0. (31) 


3B. J. Thompson, “Oscillation in tuned radio frequency 
amplifiers,”’ Proc. I1.R.E., 19, 421-437 (1931). 

4]. R. Nelson, ‘Discussion of oscillation in tuned radio 
frequency amplifiers,’’ Proc. I1.R.E. 19, 1281 (1931). 
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The second factor may be separated into its 
real and imaginary parts giving 


(QF)?+3(RwC)(QF)—1=0 (Real part), (32) 
and 
2(QOF) +3(RwC) 
+2gnR(RwC)=0. (Imaginary part), (33) 


which when solved simultaneously produces 


QF = +[(2gnR+3)/(2gmR—3)}', (34) 


and 
+2 


C= 
Rol 4(gmR)?—9}! 





(35) 


Rue 


QnR- 0 +2 
PaO \ 
mR: 2 1a! 
QmR-8 ‘ 





nw 
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t 
' 
' 
i 
’ 





Fic. 3. Oscillation conditions for a single-tuned amplifier. 
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It is necessary again to choose the negative sign 
of (QF) so that a positive value of C will result. 
As before, w may be determined from (34), thus: 


me) 


w*— wo? = a 
22mR— 


Q 


or 





2gmR+3 2gmR+378 
“rll sl Teal} 
2gnR—3 2gmR—3 
wherein the positive sign of the radical has again 


been chosen to correspond with physical reality. 
Substitution of (36) in (35) then gives 


40 





~ Reol E( (2gmR+ 3)? +40? (421° 


Again it is seen that there is a critical value of 
ZmR below which it will be impossible to make the 
system oscillate. The critical value in this case is 
gnR=%. For high stage gains and high Q’s, C 
approaches 1/wog,R?, which is the approximate 
value found by previous investigators. The 
inherently stable region is illustrated in the curves 
of Fig. 4. Here the asymptote to the hyperbola 
corresponds to the straight line with g,,R= 3 so 
that the stable region lies between g,,.R=0 and 
&mR = 3. 

The other part of the determining Eq. (31), 
(Z:+Z2)=0, produces no physically realizable 
conditions in this case and may, consequently, be 
neglected. 


Three Stages 


The definitive equation in this case is, from 
Table I, 


Z°+62722+3gmZ 7°22? + 102122? + 8gmZiZ2' 


+ gm?Z1Z2'+4Z2' + 3gmZ2*=0, (38) 


which becomes, when separated into its real and 
imaginary parts, 


(QF)*+6RwC(QF)*+[10(RwC)* —6 (QF)? 


+2RwC[2RwC —3gnR—9](QF) 
+[1—10(RwC)?—8gnR(RwC)? 


— (gmR)?(RwC)?]=0 (Real), (39) 
and 
—4(QF)*—3RwC[6+¢,R |(QF)? 
+4[1—5(RwC)?—2g,R(RwC)? |OF 
+ (RwC){6—4(RwC)? 
+3g,R[1—(RwC)?]} =0 (Imaginary). (40) 


Because an exact simultaneous solution of these 
equations is difficult to effect, an approximate 
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— —, (37) 
R?—9) }}—(2g,R+3)} 





solution is presented. If it is assumed that (g,,R) 
is very much greater than one and (g,-R)(RwC) is 
of the order of unity, then Eqs. (39) and (40) 


become 

(gmR)?(RwC)?+6gnR(RwC)OF 
+[6(QF)*—(QF)*—1]=0, (41) 

and 


(1-O@F*)[3gnR+4Q0F]=0. (42) 


It is easy to show that the only value of (QF) 
which leads to a real simultaneous solution of 
(41) and (42) is (QF) = —1. The resulting values 
of C are then 








0.764 1.526Q 
C= = , (43) 
gmR?w wogmR?[(1+407)!—1] 
and 
5.236 10.4720 
C= (44) 


~ gmR%wo  wogmR*{(1+402)#—1] 


which are the results obtained by Thompson. 
Because of the gross approximations involved 
in obtaining Eqs. (43) and (44), no region of 
intrinsic stability was found. By another method, 
it is possible to determine the stable region. If 
(QF) is eliminated between Eqs. (38) and (39), 
by some method such as Sylvester’s, an equation 
of the twelfth degree in (RwC) can be found. 
There is a theorem in algebra that states that the 
vanishing of the coefficient of the highest power 
term determines the infinite roots of the equation. 
From the previous discussion on one and two 
stages, it is apparent that this is just the condi- 
tion that is desired. It can be shown that the 
coefficient of (RwC)!* is (3gmR—4) (gmR —4) (7 gm? R? 
+16g,R—32), from which it follows that 
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gmR=1.28 is the only root having physical 
significance. The inherently stable region of stage 
gain (gnR), therefore, lies between g,R=0 and 
&mR = 1.28. 


Four Stages 


Table I indicates that the conditions of oscilla- 
tion are to be determined from 
(21° +32Z1Z24+-8mZ12Z2+Z27)(Z7+5Z1Z2 
+3gmZ1Z2+5Z?+4gmZ2*) =0. (45) 
The first quadratic on the left-hand side of (45), 


when separated into its real and imaginary 
parts, becomes 


(QF)? +3(RwC) (QF) 


+(RwC)?—1=0 (Real part), (46) 
and 
2(QF) +3(RwC) 
+(gmnR)(RwC)=0 (Imaginary part). (47) 


A simultaneous solution of these two equations 
yields 








OQF=+ ’ (48) 
[(gmR)? —_ 5} 
and 
+2 





C= ; (49) 
Rol (gmR)?—5} 


The angular velocity is found from Eq. (48) to be 
Wo | (gmk _ 3) 

gu—i{ + 
2Q' [(gmR)?—-5} 

FE — 3)? 





(gmk)? cae 5 





4 
+40 ; (50) 


where again the positive sign of the radical 
resulting from the solution of the quadratic for w 
must be chosen to correspond with physical 
reality. To produce a physically realizable value 
of C, it is necessary to choose the positive sign, 
which in turn requires that QF be negative; 
consequently the negative sign of the first 
quantity within the brackets of (50) must be 
used. Elimination of w from (49) results in 


4Q 





C= ; 
woR { — (gmR+3) +[(gmR —3)?+40?((gmR)?—5) }} 


This value of C is then the only real value that 
produces oscillation, due to the vanishing of the 
first quadratic on the left of (45). 

An exact solution of the second quadratic on 
the left-hand side of Eq. (45) is very difficult to 
obtain so only an approximate solution will be 
attempted. The real and imaginary parts of the 
quadratic in question are 


3(QF)?+ (RwC)(10+3gnR) (QF) 


+(RwC)?(5+4gnR)—1=0, (52) 
and 
(QF)*+5(RwC)(QF)*?+[SRwC — 3 ](QF) 
—(RwC)(5—3gnR)=0. (53) 


By Sylvester’s method, it is possible to eliminate 
(QF) from these two equations and obtain 
(RwC)®(S+4gmR) (mR? + 20gmR — 25) 
+ (RwC)*(81gmtR* —432gm3.R3 — 2565 2n7R? 
— 3140g,,R — 400) + (Rw()?(108¢,,7R? 
—48¢,,R—320)—64=0. (54) 
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(51) 





The first term of this expression is negligible, 
relative to the second since it is of the order of 
(RwC)*(gnR)*, while the second is of order 
(RwC)*(gnR)*. If the gain is high, then it is 
possible to neglect all the terms but those of 
highest degree in the coefficients of (RwC)* and 
(RwC)*. Equation (54) then becomes 


(RwC)*81gn'R*+ (RwC)*108¢,,2R?-—64=0, (55) 
of which the only physically realizable solution is 
C=2/3gmR?w. (56) 


t may be shown that corresponding to the above 
solution for C, (QF) = —1, so that 


w=wol —1+[1+40"]}}}. (57) 


Elimination of w between Eqs. (56) and (57) 


produces 
2 


C= ‘ 
3wol (1 +4Q?)!— 1 ]emR? 


Because of the rough approximations made in 





(58) 
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Fic. 4. Oscillation conditions for a two-stage 
single-tuned amplifier. 


obtaining these results no stable region of gain is 
found; however, the infinite roots of (RwC) are 
defined by the vanishing of the coefficient of the 
highest degree term in Eq. (54). The only 
physically realizable value of g,R is found to be 
gmR=1.18. 


Five and Six Stages 
Because of the complexity of these expressions, 
only partial solutions will be presented. Usually 
the least value of C that will cause oscillations is 
of the greatest interest so that only these results 
will be presented. By the use of approximations 
similar to those used for three stages it can be 
shown that for five stages 
0.63 


J 


1.260 


i 
ZmR*w wogmR?[(402+1)!—1] 

where (w) is determined from the condition that 

QF=-—1. The stable limits of gain lie between 

gnR=0 and g,R=1.15, approximately. 


For six stages the corresponding value of C is 


0.586 1.1720 
Co——— 2» —_____—__—_—_., (0) 
W2mR? gmRL(40?+1)!—1 Joo 
in which w is determined from QF=-—1. The 


stable limits of gain are between g,R=0 and 
gnR=1.12, approximately. 


6. CONCLUSIONS 

The general expression for the voltage gain of 
an n-stage amplifier having identical but arbi- 
trary plate impedances driven by a generator of 
any internal impedance has been obtained by the 
methods of the calculus of finite differences. The 
voltage gains, computed by this means reduce to 
the correct form as the grid-to-plate impedances 
of the tubes increase without limit. 

An application of the general expression for the 
gain has been made in determining the criteria of 
oscillation in multi-stage, single-tuned amplifiers. 
These results check previous work for large 
values of Q, but also indicate a region of in- 
trinsically stable stage gain wherein it is im- 
possible to make the system oscillate even with a 
grid-to-plate capacitance arbitrarily large. The 
stable region of stage gain lies between zero and 
2, 1.50, 1.28, 1.18, 1.15, and 1.12 for one, two, 
three, four, five, and six stages, respectively. An 
approximation, accurate within 10 percent of 
this sequence, is (2)'/" where m is the number of 
stages. 

It is apparent from the discussion of oscillation 
conditionsjthat if the actual inter-electrode 
capacity of the vacuum tubes in a cascaded 
amplifier is smaller than the least value calcu- 
lated by the results derived above, the system 
can never oscillate because of feedback through 
the plate-to-grid capacitance and can be made 
to serve as a useful amplifier. 
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The Evaporation of Antimony 


Louis HARRIS 
Vassachusetts Institute of Technology, Cambridge, Massachusetts* 
June 19, 1946 


HE aggregate formation observed in the evaporation 

of antimony has been found to have a counterpart in 
a simultaneous decrease in electrical resistivity of thin 
antimony deposits. The evaporation procedure has been 
modified as a result of these findings to produce films 
having lower electrical resistivity and higher thermal e.m.f. 
than obtained previously. 

An unusual phenomenon was observed during studies on 
the evaporation of antimony. If the heating current was 
discontinued after the antimony was first observed to be 
condensing—as noted by an evaporation ‘‘meter’’ (see 
below )—the ‘‘meter’’ resistance continued to fall during a 
period of two to four minutes. For a number of runs (but 
not for every run) the decrease in resistance with time 
was found to follow the course of a first-order reaction. 
The phenomenon is not difficult to repeat, if the heating 
current is discontinued at a ‘‘meter’’ reading of about 
5X 10° ohms (estimated antimony thickness = 5 X 1077 cm), 
the final reading of the “‘meter”’ will vary from 2000-10,000 
ohms. (It will depend on the rate of heating before the 
heating current was discontinued.) The phenomenon has 
been observed for other types of meters than the one 
described here. It is best explained by assuming that the 
meter is recording a change to crystalline aggregates of 
antimony.! 

The meter was directly (15 cm) above the evaporation 
boat (molybdenum). Two cellulose nitrate films, supported 
on mica frames, were disposed to either side of the meter. 
By looking down into the bell jar at the back of the 
cellulose nitrate films, during the supposed aggregation 
process, one can observe the growing areas of aggregates, 
as has already been reported.2* (The growing areas may 
be observed also on the walls of the bell jar.) 

We have modified the evaporating technique for anti- 
mony on the basis of these observations. We allow four 
minutes to elapse for the “crystallization” of the first 
deposit before proceeding further with the evaporation. 
The current through the boat is then increased gradually 
over a period of twenty minutes, until the boat is “red hot” 
and all the antimony has been evaporated. A vacuum of 
10-4 mm Hg is maintained. 

The progress of the evaporation was followed with the 
“‘meter’’ mentioned above.‘ The meter consists of a glass 
strip 5 em X0.35 cm. Chromium-silver is evaporated? at 
the ends and No. 22 copper wire leads are “‘soft’’ soldered 
to the chrome-silver. Leads from terminals are brought 
out of the bell jar to a ‘“‘volt-ohmist.’’ The electrical 
resistance of the meter serves as rough measure of the 
amount of metal deposit and of the resistivity of the 
metal deposit at the end of the evaporation process. Such 
a meter has been used for hundreds of evaporations and 
can be used for other metals besides antimony. The 
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antimony is removed from the meter after the evaporation, 
by solution with concentrated hydrochloric acid or by 
rubbing with a damp cloth. 

The antimony deposits are a brilliant steel-grey and 
show a pattern of tiny crystals to the unaided eye. The 
resistivity of the antimony deposits on cellulose nitrate 
films varies from 119X10-* to 60X10-* ohm-cm, for 
antimony thicknesses from 0.48 to 2.2X10-* cm. The 
thicknesses were determined from weight measurements 
and assuming macroscopic density. The thermal e.m.f. 
versus copper of the antimony deposits varied from 33 to 
38.6 microvolt per degree C. 

If the initial rate of antimony deposition is too high, 
or if sufficient time is not allowed to elapse, then the final 
deposit does not show the crystalline pattern, though it 
will still have the brilliant steel-grey appearance. Such 
deposits have a more uniform appearance; however, the 
electrical resistivity of such films is usually greater than 
120X10-* ohm-cm and the thermal e.m.f. versus copper 
may fall to less than 30 microvolt per degree C. 

* The work described here was carried out in whole or in part under 
contract OEMsr-1147 between the Office of Scientific Research and 
Development and the Massachusetts Institute of Technology. 

1L. Marton and L. I. Schiff, J. App. Phys. 12, 759 (1941). 

2 Hans Murmann, Zeits. f. Physik 54, 755 (1929). 

3 R. Bowling Barnes and M. Czerny, Phys. Rev. 38, 338 (1931). 

4 The meter used here is a modification of one used originally in this 
laboratory by A. Stockfleth (A. Stockfleth, Massachusetts Institute of 
Technology Ph.D. Thesis, November, 1942). 


_t I am indebted to Professor John D. Strong for evaporating the 
silver-chromium terminals for the meter. 
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New “Microwave Tower” 


Federal Telephone and Radio Corporation of Newark, 
New Jersey, is building a 300-foot ‘“‘microwave tower’’ for 
Federal Telecommunication Laboratories at Nutley, New 
Jersey. It is believed that such towers will be headquarters 
in cities and towns for simultaneous use of many microwave 
radio functions, including f-m broadcasting, Pulse Time 
Modulation broadcasting, television, police radio net- 
works, communication with mobile units, radar applica- 
tions, aerial navigation, and the interception of illegal 
transmissions. At the ground-breaking ceremony Colonel 
Sosthenes Behn, President of International Telephone and 
Telegraph Corporation, stated, ‘‘Various broadcasting 
services in important communities will find it much more 
efficient to cooperate in building one great tower to be 
shared by all functions utilizing microwaves than for each 
service to build its own tower. Joint financing would make 
it possible to build a tower superior to any that could be 
afforded by individual services.” 

Federal’s microwave tower will be an experimental 
laboratory, first of its kind ever to be built. Microwave 
experiments heretofore have been conducted in tall office 
buildings and other structures not well adapted to the 


purpose. 
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Ground has been broken also for the newest addition to 
Federal’s existing Nutley laboratories, which will increase 
the present floor space of the laboratories by 65,000 square 
feet. 

Meteorology Experts Available 


The American Meteorological Society, through its 
president, Professor H. G. Houghton of Massachusetts 
Institute of Technology, wishes to call to the attention of 
American colleges the fact that a large number of carefully 
selected college students were given training in professional 
meteorology at a few selected universities in this country, 
preparatory to their employment by the Army and Navy 
Meteorological Services. A number of these men, before 
entering the meteorological field, had done a considerable 
amount of graduate work in physics and are now available 
for teaching and research in both these fields. 

In most of our colleges which offer non-professional 
courses in meteorology to meet the needs of students in 
geography, geology, and agricultural subjects, it has been 
standard practice to have these non-professional courses 
taught by members of various departments. The return 
of many of the professionally trained service meteorologists 
offers an opportunity for our colleges to obtain the services 
of instructors who are capable of doing professional work 
in a physics department and at the same time are qualified 
as experienced and trained professional meteorologists to 
present stimulating non-professional courses in meteorology 
on a somewhat higher plane than may have been possible 
in the past. 

The Society is anxious to cooperate with American 
colleges in providing information concerning these re- 
turning weather officers that might be helpful in placing 
them in suitable teaching positions. To that end, files have 
been kept which contain detailed information on the 
education and practical experience of these men. Inquiries 
may be directed to the Executive Secretary, American 
Meteorological Society, Five Joy Street, Boston 8, Massa- 
chusetts. 


Books for Europe 


The American Book Center for War Devastated Li- 
braries, Inc., is collecting and shipping abroad scholarly 
books and periodicals which will be useful in research and 
necessary in the physical, economic, social, and industrial 
rehabilitation and reconstruction of Europe and the Far 
East. The Center cannot purchase books and periodicals; 
it must depend upon gifts from individuals, institutions, 
‘and organizations. 

What is needed: Shipping facilities are precious and 
demand that all materials be carefully selected. Emphasis 
is placed upon publications issued during the past decade, 
scholarly books which are important contributions to their 
fields, periodicals (even incomplete volumes) of significance, 
fiction and non-fiction of distinction. All subjects—history, 
the social sciences, music, fine arts, literature, and especially 
the sciences and technologies—are wanted. 

What is not needed: Textbooks, out-dated monographs, 
recreational reading, books for children and young people, 
light fiction, materials of purely local interest, popular 
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magazines, popular non-fiction of little enduring signifi- 
cance, etc. 

How to ship: All shipments should be sent prepaid via 
the cheapest means of transportation to The American 
Book Center, care of The Library of Congress, Washington 
25, D. C. The Center cannot accept material which is 
sent collect. When possible, periodicals should be tied 
together by volume. It will be helpful if missing issues are 
noted on incomplete volumes. 


New Section at Midwest Research Institute 


The Midwest Research Institute, Kansas City, Missouri, 
recently announced the formation of its Engineering 
Mechanics Section. Previously, work in this field was done 
under the supervision of the Applied Physics Section. 
The new section is equipped with a large physical testing 
laboratory, machine shop, and electromechanical labora- 
tory. Problems within the scope of the group include 
structural and mechanical design, vibrations, fluid flow, 
electron application in engineering, and servo-mechanisms. 

Martin Goland, formerly of the Structures Department 
of the Curtiss-Wright Research Laboratory, Buffalo, New 
York, has been appointed chairman of the new section. 
Clayton O. Dohrenwend, formerly chairman of the Me- 
chanics Departments at the Illinois Institute of Technology 
and the Armour Research Foundation, will serve as re- 
search consultant and senior scientist. 


Food Packages for Europe 


The Cooperative for American Remittances to Europe, 
Inc. (known as CARE) has announced a new service for 
sending food packages to Europe. CARE is a non-profit 
corporation organized by 24 relief agencies, with the 
support and approval of the President’s War Relief Control 
Board. 

CARE accepts orders for individuals, groups, or organ- 
izations in Austria, Czechoslovakia, Finland, France, Italy, 
Netherlands, Norway, and Poland, and recommends the 
purchase of a 10-in-1 package of Army rations for $15 
each, with delivery overseas guaranteed. Orders when 
received here are rushed to Europe and filled overseas from 
a stockpile of CARE food packages already over there. 
When the food is delivered, CARE’s overseas representa- 
tive will get a signed receipt which will be sent to the 
donor. If a package cannot be delivered within a reasonable 
length of time, CARE will refund the full purchase price 
to the donor. A typical CARE package contains such 
products as Borden’s butter, Nescafe, Swift’s canned ham, 
Del Monte corn, Welch’s jam, Hormel’s sliced bacon, etc. 

CARE’s order blanks are available at banks and are 
also distributed by relief, fraternal, professional, and 
social organizations, but if you have difficulty in getting 
one, write to CARE, 50 Broad Street, New York 4, 
New York. 


Annual Convention of I.E.S. 


The Illuminating Engineering Society goes to Quebec 
September 18-20 for its Victory Convention, sharing with 
its Canadian membership this annual event which during 
the war years was limited to technical conferences. 
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